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PHYSICAL OF DRAIN 


te) 
SYNOPSIS" 


as slope of the principal drainage channel, sor of the princi- 
pal drainage channel, shape of the drainage-basin (basin shape), and general b 
Slope of of the land in the basin (land slope), are reviewed and defined in terms 


. the standard ; area- -distance curves and the e modified area- -elevation curves are 


the engineer is lucky if there is astream- -gaging 

in the immediate vicinity upstream or downstream of the point in which he - 
interested. Thus in most cases, he is required to estimate stream ! flow by com- 

paring the drainage - -basin above the point of interest with another basin 
similar physical characteristics. In the case of a highway, the hydraulic engi-- 
neer generally computes discharge corresponding to a particular frequency by J 
a for that frequency and entering with an ‘argument in 


_ ciety of Civil Engineers, Vol. 86, No. HY 3, March, ee a I 
Head, Hydr. Dept., Howard, Needles, Tammen & New York, 
ay: Hydr. Engr., Howard, Needles, Tammen & Bergendoff, New York, N. Y. 
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‘ 
OF recent Usape. Stahdare cedu p Ang the ar qs pons 
and area-elevation curves is suggested. Also discussed isa me on 
ifi - ion curve. Numerical pro- 
veloped for computation of a modified area-elevation 
— 
— 

| 
= = a written request must be filed with the Executive Secretary, ASCE. This paper is part ey. te 
. — 
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_ have been developed by y correlating discharge against - physical ‘characteristics 
of many similar drainage-basins. In the case of a dam, the unit-hydrograph + 
method is used. Coordinates of the peak value of the unit hydrographare taken | 
4 froma graph in which the peak (maximum y- coordinate) and lag (corresponding z 
x-coordinate) of unit hydrographs for basins havi ing similar physical charac- 
i teristics have beencorrelated with the physical characteristics of those basins. 
— Many recent papers have been published using such correlations so that a 
knowledge of basin characteristics | has become an essential for any — 7 
== various physical characteristics of drainage-basins have t been defined 
‘in ymany different publications in in the past by many different authors. 3, 4,5, This 


to such publications in many recent papers quite useless. iQue — 
In an attempt to remedy this situation, the various drainage - -basin charac-_ 


THE DRAINAGE ‘BASIN 


is generally in square miles (or for small areas). 
ally, sub-surface drainage - basins (ground- -water flow) encompass different 


areas. _ However, this complication is very rare —_ generally not of great 


The length of the principal drainage channel is measured 


length is s normally measured a as s short chords | on n United States Geological Sur- 

_ vey (USGS) quadrangle maps, which along with the projection of the principal 
channel to the basin divide, may cause some minor discrepancies in the length 
_when measured by different individuals. However, such discrepancies are 


enerally insi nificant in their effect on hy drolo ic com utations. ety 


= “Drainage Basin Characteristics,” by R.E. Horton, Transactions Amer. - Geophysi i- 
ay 5 “Synthetic Unit Graphs,” by F. F. Snyder, Transactions Amer. Geophysical io 
Vol, 19, Part 1, 1938, pp. 447- 
_ 6 «Unit-Hydrograph Lag and Peak Flow Related to Basin Characteristics,” by A. B. 
"Taylor and E. Schwarz Union, Vol. 33, 1952, PP. 


| 
> 
3 
L 
a= 4 several different ways. Also, many of these original publications are no long- — — 
— 
— 
| 
‘The area O € drainage-baSin 18 Its most Importan physical character- 
istic. The area is actually the horizontal projection of the land surface 
= j which run off into the surface channels, above the point of interest, occurs. — It eee 
is ordinarily measured on topographic maps by location of the divide, or sepa- 
reel rating ridge, which determines whether surface run off will flow toward the as 
— “ee 
— — 
— 
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DRAINAGE BASINS 


_ The slope of the principal drainage channel is probably the most significant _ 
physical characteristic after the area of the drainage-basin. ‘The simplest 
method for expressing the slope of the principal drainage channel is to divide _ 

, ~~ length of the channel by the difference in elevation between its upper —— 
— ends (Definition 1, Fig. 1). Generally, the length of the channel is mea- 
- sured beyond the upper end of the clearly discernible stream channel to the | ni 
_ drainage divide and then divided by the difference in elevation between this — 
a point on the ridge line and the point of interest on the channel. The USGS has _ 
used a slope parameter instudies inthe New England states in which the length — 
_ of that part of the stream between points 85% and 10% of the total distance above 
_ the point of interest is divided by the difference in elevation at these places. - 


—_, are of th the e opinion that the most st upstream part of the slope, in the seep 


ELEVATION - DISTANCE CURVE> ¢ 


SLOPE DEFINITION 3 


Draw line so Area > 
above channel profile 
Area below 


ELEVATION 


Principal 
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, —SLOPE OF THE PRINCIPAL DRAINAGE 


nished by the headwater area, and also that the most downstream part of the 
~ slope of the principal channel may be so flat that ‘it might | not indicate the true 
portion of the flow which the lower | portion of the drainage area actually « does’ 
A slope of the principal drainage channel can 1also be determined by dividing 
_ twice the mean height of the channel profile above _— 7 of rie: by the 
ah _ The mean height of the profile is obtained by integration of the elevation- 
distance curve (channel profile) and dividing by the length of the channel (Defi- 


nition 3, Fig. 1). . The integration can be accomplished by “ “planimetering” the 
area below the elevation-distance curve. This slope is also _— by the = 


= 
7 
— 
— 
— 
= 
— 
” by M. S. Benson, Proceed- 
“Channel-Slope Factor in Flood Freq ASCE Vol. 85, No. HY 4, April, 1959, pp. 
al of the Hydr. Div., ASCE, Vol. 85, 
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of the line drawn through the « origin on stream- curve such that the 


; areaunder it equals the area under the profile curve. The mean _ channel height — 
is at a pointon the above-cited line at one-half the lengthof the principal drain- | 
age channel measured from either end. Rather than computing mean channel | 
height as such for use in determining channel slope, the mean basin height (to _ 
be defined, subsequently, under the heading “Land Slope” )is often used and the é 
resultant slope is called the mean basin slope. - This effectively introduces a 
weighting factor to account for slope of the basin aswell as slope of the stream 


ames channel is the slope of an equivalent stream having the same travel 


and same length (Definition 2, Fig. 2. This slope is computed 


where p equals the anaes r of equal nin. ‘into es the stream has been 


the summation of the( values for each reach. 40 Ea. ¢ is on the as— 
sumption that velocity is proportional to the aquere root of the slope. This 


_method of determining slope was defined by A A. B. Taylor and H. E. Schwarz, 


M. ASCE and is used by the Corps of Engineers. ERA E 


The latter two methods of slope determination give values that are 


mately equal an — which may for small (Fig. 1 
BASIN 


The shape pe of a -basin can ed in terms of various easily, 
> defined and measured factors. One measure of basin shape is given by the 
quotient of the drainage area divided by the length of the major watercourse, 
_ - A/L. Another measure of basin shape is given by the length of the major water- 
course squared divided by the basin area L2/A. 
; The most commonly used basin shape factors are Lea and the product L | 
. times Leg to some power, , usually (L Leg)0.30, where L is the length along the 
ar longest watercourse and Lea is often defined as the distance along the main 
aca | Gralange channel from the point of interest to a point opposite the computed © 
center of | gravity (centroid) of the drainage 2 area.9; 6 A more accurate definition 
of Lea is given subsequently. The term lea | is computed by integration of the yo 
area-distance curve and dividing by the drainage area. 
- Toconstruct the area-distance curve (Fig. 2) the channel length is first sub tae 
divided into separate reaches. These reaches may be of equal length along the 
watercourse, such as 1-mile reaches or equal parts of the total length. Divid- P 
ing points along the stream couldalso be made at points where contours cross. yo 
A preferable method, especially for larger rivers, is to divide the watercourse 
: length just above and | below e; each major tributary. en _ The areas in this c case are 
- the contributary drainage areas at ‘these | points . After: the channel is divided 
_ into reaches, the sub-drainage basins contributing to the stream between the 
limits of each reach are determined 


ae 
ae... 8 Unit Hydrograph Compilations, Project cw 153, Washington | District, Corps of 
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E BASINS a 


, = A curve of channel distance upstream from the the point of interest (or per- 
centage of channel distance) v versus area (or p percentage of area) contributing 
Within that distance is then drawn as shown in Fig. 2. The distance to the cen-_ 
troid of the area is found by dividing the area above the curve ve by the drainage © 

- area (or percentage of the drainage area) which is the ‘maximum ordinate . AS 
an alternative, Leg may also be found by integrating the area-distance curve | 
and dividing by the drainage area. The integration is accomplished by  “plani- — 

metering” the area above the area- -distance curve. 
-¥% A suggested method for plotting the area-distance curve and a sample ane 
of Leg based on this curve will be presented subsequently. 
Physical Definition of | a-— Previously, it was stated that the centroid of 


_ the area above the area- distance curve determines Seats In other words, —- 


is the length measured up the stream channel from the base of the drainage 
_ area toa point corresponding to the the centroid of of the the area above the area- ona fll 


Orainage Area 
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ad 
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DRAINAG 


DISTANCE uP PRINCIPAL CHANNEL 
«FIG, 2. —BASIN SHAPE FACTOR, 


curve. all of the e various methods of approach ‘that utilize the 
function derive it in the same general manner, some writers have added physi- 
is definitions of lea which are not rigidly correct while others leave the 
physical interpretation of Leg purposely vague. For example, Leg has been 
defined as “the distance along the stream channel to a point opposite the cen- 
troido of the ‘drainage area.”6 This ¢ definition is actually incorrect since the lo 
a cation of the centroid of an area depends solely on its shape, and L<a, derived . 
> from the area-distance curve, obviously depends on other factors as well - On 
hand, Lea is sometimes vaguely definedas “the distance to the center 
Physically speaking, just what is ? Fig. 3 shows sa drainage 


i area. Axis x-x is drawn perpendicular to the major axis of the area through 

OM “Applied Hydrology,” by R. K. Linsley, M. A. Kohler, and J. L. H. — 


| 
— 
| 
ROM BASE 
— 
— 
| 
— 
il: 
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= in which AA is an elemental area, y] is ‘its menses from the axis x-x, and A. 


a is the total area. The e term Lea can be expressed by a similar formula 


7 . in Eq. 3 y2 is the distance | measured along the stream to the point on the ~~ e 


Lore 
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De Defining a characteristic which adequately describes land slope f has proved 7 
‘difficult. However, this characteristic is a major factor in indicating the travel 
time of overland flow, so several methods for determining land slope have been — 


Most of the methods generally used for defining land slope take into account — 
se general slope of the entire basin including the tributary stream network :. 
contributing to the major watercourse, but excluding the major watercourse. = ; 
Average land slope | has been defined as equal to the plotted contour interval 
times the total length along all contours in the drainage-basin divided by the © =a 
. drainage area. However, measuring the length of contours proved a tedious 7 


- and time-consuming operation and hence the following short-cut methods * 


— 
7 
4 


In the intersection- line method, a grid of siti spaced parallel and 
perpendicular lines is laid over a contour map of the drainage area. The num-— 
ber of contours crossing each s subdivision of the oid within the areais counted, _ 


‘The land slope in either grid direction is then es = where N is the total : 


_ numberof contour crossings for all lines in one direction, AZ denotes the con- 
4 tour interval, and L is the total length of grid lines in one direction within the ri 


where N’ and L' are the sum in directions of ‘the values 
defined as N and L, and @ is the average angle of intersection between the con 

= and grid lines. yi Often a value of 1.57 is used for sec 6 as it is the av-- 
erage secant of angles from 0° to 90°. The apy of ee uses this method — 
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-AREA - ELEVATION CURVE 


DRAINAGE AREA A 
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_ Another method of defining land slope, or more accurately, basin slope, util-— 
izes an area-elevation curve in which the drainage area above a certain ele- 
vation is plotted against that elevation. Area-elevation data can be assembled 

_by planimetering the area enclosed by each contour and the basin divide and 

plotted as shown on Fig. 4. 


er The mean basin height can be obtained by integration of the area- elevation — 


curve and dividing the result by the drainage area. It can also be obtained by 
the method shown in Fig. 4 in which a line is drawn such that AL= ne + Ag. 


the drainage a area “indicates the ‘mean basin height. . Mean basin slope is ag 
defined as two times the mean basin height divided by the length of the channel, 


— 
7 
ihe average lana siope is wen given DY some relation between these two 
res 
— 
— — 
— — 
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Re Another curve can be plotted conveniently by using the same data derived 
for producing an area- -distance curve. This is called a modified area- ~elevation 
curve. Rather thancutting the drainage- -basin into slices between adjacent con- 
tours, it can be divided into sub-drainage basins contributing to equal reaches 
along the length of the principal drainage channel to correspond with the method 
used in computing the area-distance curve. ‘The preferable method, especially 
for larger rivers, is to break the channel at the tributaries just above the point 7 - 
_ where they jointhe main channelas described under the heading “ Basin Shape.” 

In In this case, area no longer refers to area lying : above a certain elevation but : 
= rather to the area contributing 1 runoff to the major stream above the Polat at 
_ which the stream bed is at a certain elevation. Thus the mean basin height is 
obtained actually weighted by distribution of drainage area long the stream 
file. pI Mean elevation and mean basin slope are correspondingly weighted. - 

The foregoing method saves considerable computation time over the 
method for computing an area-elevation curve by utilizing one set of area data 
_ as the basis for both curves . Inthe opinion of the writers this results ina 
‘more ‘meaningful mean basin height and mean basin slope. it is undoubtedly 

not applicable in studies of snow hydrology where the temperature- may ate . 
curve makes true land elevations more significant. _ Further study may also 


SUGGESTED PROCEDURE FOR DEVELOPING CURVES 
In the plotting of the area- curve to compute mean ‘elevation and 
_ mean basin height and the area-distance curve to compute ca it makes no ; 
- difference whether zero channel distance is taken at the base of the drainage 
area or at the drainage divide. Also it does not really matter whether the 
. area above a certain point on the channel or the drainage area below A 
_ porn vary in any one of several ways. However, in order to eliminate any : ; 
possible confusion a suggested standard procedure similar to the method used 
by the Corps of Engineers8 is outlined herewith: 


P a Channel distance will always be measured up the channel from the base and 


a 


is 


ill be expressed in percentage of total length. Elevation will always be mea- | 
_ sured above the stream-bed elevation at the base of the drainage- basin and will © 
be expressed as the percentage of the total elevation drop from one end of ol] 
_ Stream to the other. Drainage area will also be expressed as a percentage of 
ee the total. It will be recalled that in the development of the area-distance re-_ 
__ latlonship the term area referred to is the area of the drainage -| -basin contri 
°< ing below the point along the stream in question. On the other hand the a 


elevation relationship required the consideration of the contributary area above 


the point in question. This distinction need cause no > confusion in actual prac- | 


oo 7 The fact thatall scales’ go from 0% to 100% alstiie the engineer to o plot both | 
_ curves on a square piece of graph paper. Since both curves appear roughly as — 
diagonals on such a graph, the scales are oriented so that they will cross only 
- once, thus” eliminating any danger of confusion (Fig. 5). . The value | of Lea is 
4 then found by measuring the area above the area-distance curve in square inch- - 
and to ‘square percentage, for for verti- 


| 
4 
a 
aoe 
: 
— 
— 
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DRAINAGE BASINS 
cal scales of 1 in. = 10%, the conversion factor is 1 sq in. = (100%)2. “This | 


7 figure divided by 100% is the Leg expressedas a a percentage « of the total length. 7 
Tofind mean basin height (mean elevation), the area above the area-elevation | 
curve is determined, converted to ‘Square percentage | and divided by 100%. The e, 


Mean basin ‘slope is : then mean n basin height divided by one-half the total chan- 2 
nel length, both values having been converted from percentage to actual lengths. _ 
= The modified area-elevation curve is plottedin a similar rfashion to that d de- . 
scribed for the standard area-elevation relationship. _ 
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DISTANCE (rom gouge) on 
FIG. 5. HYPOTHETICAL DRAINAGE BASIN 
Example.—The necessary points for plotting a modified area- elevationcurve 


and a an — ~distance curve for | a 1 hypothetical drainage- ‘basin are shown 


a 
c 
a 
w 
oO 
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“Note that distance i is measured in an upstream direction ‘and that area and ele- 
vation are "measured below the point on the con The modified area~ 


— 
— 
"3 | 
— 


eine above, which is determined - subtracting percentage area below ‘iam 


100%. The modified area-elevation curve can directly using the 
origin i in the upper-right corner. 


TABLE. —HYPOTHETICAL DRAINAGE BASIN (pRADUGE AREA = = 8.10 


"CHANNEL LENGTH TO BASIN DIVIDE = 5.42 MILES) 


"Elevation, 


Pees 


0.1270 
0.0878 
=: 0.00735 | o.osss | 4 


( = = 0.0061 ,00615 ft per ft & 


127.5 


A ‘sample computation of the weighted values of mean basin height, n mean 


=. elevation and land slope (mean basin slope) based on the data given in Table 1 


modified area-elevation curve = (38105)? 
Mean Basin Height = (0.381) x (221 - 21) = 76. 
Mean Elevation = 21 + 76.2 = 97. 2 ft (m.s.1.) 
‘Land Slope (Mean Basin Slope) = = 310 = 0. 00532 ft per ft 
28. 1 ft mile 


= 

| ae | | se | ot 

68 9 | 2860 | 0.00315 | 0.0561 — 
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DRAINAGE 


_ Area above area-distance curve = = = (433 2 


Lea = (0.433) x (5.42) = 2.34 miles 


Deteveinedton of the Channel Slope.—In Table 2 stream slope for the same 
: hypothetical drainage basin ¢ considered in Table 1 is computed by the method | 
_ previously reviewed in which the slope of the principal drainage channel is ex- a 
_ pressed as rr a of an equivalent stream having the same travel time and 
Same length.§ 8 Pertinent information from Table 1 ts thus repeatedin Table 2 
_ SUMMARY 
hes been presented along with a suggested standard procedure of plotting the = - 
area-distance curve, from which Lea can be computed, and the: area- 


is 


The computation of Le. is as follows: 
a 
Anadditional characteristic, developed from the modified area-elevation curve, 

— 
— 
— 
— 
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By David 


A . theoretical study is made of the development of velocity with time when 


i a a valve at the end of a drain pipe is ‘suddenly opened. The - effect of friction, 


% tained by considering the fluid inertia, is compared to the velocity obtained - 
the more _well- -known method which disregards” the fluid inertia and assumes 


“the liquid from the tank to ascertain that it is sufficiently high to give rapid _ 
: a reactor shutdown. The mechanics of rapid liquid-level change or of rapid tank 
—, falls in the realm of non-steady flow. The moderator or reflector 


% a written request must be filed with the Executive Secretary, ASCE. This paper is part 
4 of the copyrighted Journal of the Hydraulics Division, Proceedings of the ean — 
ciety of Civil Engineers, Vol. 86, No. HY 3, March, 1960. a Se ae a 
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dealing with the shutdown Of a reactor by the rapid draining Of moderator or 
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peak velocity is the same ¢ order of magnitude of time, ‘from a nucl ear point of of 
view, in which the shutdown must be completed. cieamaeaal —_ 


The more well- known type of tank-draining problem is the one in _—- -— 


but rather one of slowly varying ‘steady flow. In flow under decreasing head 
the inertia forces are neglected, whereas in non-steady flow the presence a 
a inertia forces is implied. ‘The flow relationships that are derived without the : 
inclusion of inertia forces are designated as “terminal flow” relationships in 
this paper, since they — an approximation of the terminal phase of aadl 


us vertical 


Horizontal length 
of drain pipe 


related the — of a conservative mass 
nd energy system has already been studied. The present problem deals 
q with a non-conservative mass and energy system, in which mass and energy 
in the form of draining liquid is continuously being lost from the system. we 
The oe consists of a large pressurized tank from which liquid drains 
= eee drain pipe of smaller diameter cae a vertical | anda horizontal — 


— 
“Hydrodynamics. of Liquid Poison Scram System, Burgreen, Nuclear 
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TANK ‘DRAINING 
“tank and drain | drain pipe summation of forces ont the in the tank 


It is assumed in 1 1 that friction forces retarding the flow in 


in which xX represents the displacement of a fluid particle in the drain 1 pipe | 


¥ 
is most the continuity relatio relationship AY y= = A2 x = The 


In Eq. 3, + F represents the of lost 
through pipe friction and would include entrance, exit, bend, and ‘valve losses a, os 
addition to the straight pipe losses. At the. junction of the pipe | and tank 


"provided the length of the drain pipe is greater than its diameter. —_ 


The use of the > continuity relationship, Ay ie Ag? Eqs. 4 3, 
Ag a sf. Ag 


5 is expressed in terms of the ‘displacement of a in the drain 
_ tube. The flow could be expressed just as readily in terms of the drop of level _ 
it int the > tank, y- In In 1 deriving Ea. 5, Pop the pressure the e end of 


‘At the 


ama ACCELERATION 


a 
= 
iii 
4) 
} 
* 
2 


where R = and n Po} pressurization in the tank i is thus ex- 

_ pressed as the number, n, of initial tank heads. The ‘initial acceleration | of the 

: liquid in the which is the same as the of the sur- 


t=0 
is clear from Eq. 
the liquid in a tank, to which a drain pipe is attached, can never be as high as 
| oa R <1. ts In order to obtain an initial acceleration of g, the require 


AREA oF A VERY SM. DRAIN PIPE 


L 


(n+1) H+ Ly] = 0° 


H + Bhs 


(2 at? L+RH 


and represents flow under constant head with inertia | effects. It is thus intica- a 
tive of the initial flow behavior only, since the approximation of Eq. 9 is not 
= for large x. The first integrals of Eq. 9 with respect to time and dis- ’ 


placement are, respectively, 


at | 


and the displacement may be wha explicitly as a function of time me by com- 


{2 g[ (n+) H+Ly ](1+F)} (11) 
‘If there were a ‘horizontal: drain pipe only and no no pressurization or ne - 

1 
Eq. 12 has been derived3 from the Bernoulli _ equation for non- steady flow. — 


From Eq. 10b— it is seen that for small values: of x/L the v velocity in the 
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4 a constant initial acceleration. When the drain pipe area is extremely small, - 
-_ the change in head is very slow and the starting velocity of the terminal flow 
is not sensitive to fluid movement in the — drain pipe, and can be ex- 


problem in which the is governed by the 
a ae head and the inertia of the fluid to a change in velocity is disre- 


es ‘ very close approximation of the flow for the case when both the area vatio 2 
7 is small and the length o of ‘drain pipe is small. _ Under zero inertia conditions, — 
the ‘solution of aid 5 in terms c of the - y-coordinate, or liquid level in the tank, is 


dt 


The integral of Eq. +15 gives the ‘required for a in liquid 
[ (n+1)H+L, - Yo) nis (n+1 )H+ 
"where Yo } is liquid level at =0 


5 


Although Eq. 5 is a non-linear differential equation with variable ‘coeffi- 
cients, the first Sage, in terms of the tank level y is seatnny obtained. It 


(n+1)H+L,+ 


‘effect these parameters will be ‘examined separately. first the 
_—— when the drain pipe lengths Ly = Lg = H, when there ov no — 
gation, and the friction are Eq. 17 becomes oat 


_— 18 is : plotted in Fig- 2 for several ratios of drain pipe area to tank area. ™ 


When this ratio is unity, the configuration PARRA a ee elbow of uniform a: 


i 
iS 
— 
— 
4 
— 
— 
— 
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= 


diameter whose is “the horizontal length. Eq. 18 then 


Fig. aa free fall is included for 
hows the velocity of a free- falling body, from a H, at ele- 


EL IN "TANK. 


Iti is a special case of Eq. . 20, or the R- = when Lo, the hori- 


zontal length of pipe, is equal to zero. The curves in Fig. 2 show the rate of 
fall of the liquid surface in the tank to decrease with decreasing area ratios. 


er ‘The range in which the flow is governed arlene? by inertia forces also be Be 


comes smaller with decreasing area ratio. This may be observed through the | 


increasingly early blending of the true flow curves with the terminal flow 


curves 


‘curves: as the area ratio > Saree ye al The dotted terminal flow curves 


— friction and pres- “a 
— gurization, for the disc length Lo 
total length L + H and horizon 

if 

— 09 1.0 il 


-displacement curve shows that ‘the inertia’ forces 

mine the character of the flow. - For~ large area ratios the velocity- 
> ‘displacement curves resemble the free-fall curve and any approximation with 
flow would not be valid. that the true velocity-displacement 


curves indicate, properly, that the starts at zero and builds up, 


(0.2 03 04 05 06 0.7 08 0.9 


(hing FIG, 3. —VARIATION OF FREE SURFACE VELOCITY WITH LIQUID 


_ LEVEL IN TANK; PRESSURIZATION: ‘n=10; F 


LN TANI =0; L}=L,=H. 
‘The effect of pressurization is shown in Fig. 3. The pressure on surface 
a é the liquid is assumed to be constant and is arbitrarily taken as ten tank 
_ heads; that is, n = = 10. These curves show that a rate of drop in liquid level, 
in excess of the free- -fall rate is easily obtained. It will be observed that for 


_ the selected pressurization, an area ratio of (1/4 promeces a rate of = in 


ig 
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March, 1960 
level approximating free fall. The terminal flow flow approximation is is now 
rather poor even for small area ratios. un 
_---' The result of taking friction into account is demonstrated in Fig. 4. ‘The 
total friction loss in the system is arbitrarily taken at three velocity heads — 
4 (drain-pipe velocity). 
_ and 3 it will be noted that the \ velocities, for “corresponding az area ratios, : are 
smaller and that terminal flow is established sooner. Although terminal flow 4 
is now a better approximation to true flow, ‘it is obviously not a satisfactory | 
approximation, for moderate area ratios, of the early stage flow. | en nid 
_ The effect of drain pipe length on the development of flow is well demon- 
strated by a configuration having a a area ratio atio typical oi 
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G. 4, —VARIATION OF FREE SURFACE ‘VELOCITY | with 
LIQUID LEVEL IN TANK, Ly = Lo =H; .n=0; F=3. 


wate water | with drain pipe attachment. An area ratio of 0.01 is selected, and 


the velocity- -displacement relationship is plotted in Fig. 5 wal lengths of drain 
ar pipe of zero, H, and 5H. The curves show, in common, that the inertia 
7 _ transient is completed and terminal flow started when the maximum rate of a 
_ The curve in Fig. 5a represents the ‘development of flow when a very long — 
; drain pipe, equal to 5 H, is attached to the tank. The peak velocity is not at-_ 
; tained | until about half the tank has been drained. This indicates that a long 
g drain pipe produces a long inertia transient, and that computations of dis- h 
charge based on the head lead to incorrect results - — even 


— 
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_ area ratio is small. When the drain pipe is of moderate length, equal to H, 
‘ then the peak velocity is reached when one-tenth the tank has been drained. | 
_ The fundamental differences between inertia flow and terminal flow are > 
brought out by the curves in Fig. 6 in which the fluid acceleration is plotted 
_ against the drop in liquid level. Since terminal flow is flow without oniere- . 
tion, the deviation from terminal flow is measured by the magnitude of the 
acceleration. . The parameters are the same as those used in the velocity- _ 
displacement curves of Fig. 2; namely, unpressurized frictionless flow with 
‘Le = H. The ‘uppermost a } horizontal line representing 


4 


0.01 


fall, or the equation = =1. An examination of the curves will 


thereafter continuously with increasing displacement. As the 
area ratio becomes small the curves show a range of decelerating flow which 
7 _ corresponds to the region of decreasing Terminal flow is repre-- 
sented by the abscissa, (1/g) (d2y/dt2)=0. The general trend of the curves 
is toward zero acceleration as the area ane approaches zero. - This general 
trend will prevail also for configurations other than Lj = L2 =H, but with 
e 


longer drain pipes will hoes toward zero acceleration at a slower — he 
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= obtain time-displacement curves it is necessary to integrate Eq. 17, 
Analytically, it appears to be possible to perform the integration only for spe- 7 
cific values of F and R. In general the time-displacement curve would be ob- 
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ths the expression for velocity when the acceleration is constant. This ap- s 


proximation may be used to start for accelera- 
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FIG, 7. _CONFIGURATION OF DRAIN SYSTEM 


d pipes, itis expedient to 1 use an annular weir micas as shown in Fig. 4. ¥ In ' 
this arrangement the length of the drain pipe L = Ly + L2 is assumed very 
small in comparison to g/ApE. ‘The solution of Eq. 5, disregarding pres- 7 


1-2- 
Eq. 25 is plotted i in in Fig. 8 ‘or ratios of pare pipe area to tank area of 1/2, 


1/4, = corresponding terminal- flow curves ‘are also shown. Note that 
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drain time is 
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as 


This is shown in Fig. 9. a 


flo 
2 


“FIG. 8 


— ows 
application, would give a 
now, in the = rue OW, an ora prac 1ca p drain a tank. It is, 
ee = proximation 0 “8 f the time required to completely « the starting flow, [ii 
© hen the interest is 
— 
| 
— “| 
a | 
| &§ 
ome 03 04 OS 06 0.7 08 09 |. 
— 
G 


“that about three tenths of the tank has been drained the time - 

Cc is practically t that of free-fall. The divergence becomes more 

significant as the draining of the tank nears completion. complete the 

- draining, 25% more time is required than the time it takes for a body to fall — 

- _ freely through this distance. - Terminal flow is not at all representative of the 

true state of affairs. The fact _ the free-fall and terminal flow oa meet 
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At the uppermost point of the drain pipe, near the junction with the 
the pressure in the liquid may be quite low at some stage of the flow. If the 


at this” ‘Point than vapor of the liquid, then flow 
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Por is vapor pressure of the ‘hen (Pog/ 
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pepresent of a of liquid that can be supported when separa- 
is a about to occur. Separation will not occur as long as (Per/W) < 


“8 dt2  2¢ \ at/ 
and separation will occur occur when 
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RAINING 
or to occur, rite wou 


(32), 
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At any given or level displacement there will be 
so ciated acceleration and velocity. When the acceleration and velocity tail t to ; 


place. When ‘takes place | the draining of the tank is no ) longer 
governed by Eq 
_ Before the pressure at any point in the pipe reaches the absolute vapor — 
pressure of the fluid and produces flow separation, it will pass through a pres- — 
sure range less than atmospheric. There would then be a tendency for air to ¥ 
back = through the pipes. . However, as in the case of a siphon, as long as_ 
a substantial flow velocity in the pipe al will not 


— the flow and the » pipe will flow full. 7A 


| 
The n of the foregoing of tank draining the 


4 ‘When the ra ratio of drain pipe area to tank ik erdie- sectional area is email 
ei ~ and the length of drain pipe is short, then the conventional manner of retell - 
ing the rate of draining, which disregards the inertia of the fluid, may be used ay 
_ 2. When the area of the drain pipe isa rR fraction of the tank ; area, nll 
rate of draining is given by Eq. 17. The neglect of the fluid inertia, in this 
3. When the length ¢ of the drain pipe is a large fraction o: of, or r greater than 


the initial nang of fluid in the tank, it is necessary to use ~ 17 to or 


Since either drain have a a strong influ-_ 
ence on the flow, a combination of moderate area ratio and drain pipe length | 
also requires Eq. for the proper computation of the rate of draining. 
§. As one may expect, ‘the effect of drain- -pipe friction is to decrease the 
: “g rate of draining while the effect of tank pressurization is to increase the rate 
= “of draining. For expediency the form of the friction term was taken | as the 
3 om the flow is actually 1 non- -steady. 


"APPENDIX. DERIVATION OF MOMENTUM EQUATION ‘BY 


when dx/dt = 0, in order ould be necessary that 
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‘he potential energy, using asa datum is 


for the ‘drain p pipe. 
lost from the system, in the form of is = 
aw 
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2 Agta) 
given by Eqs. 33 through: 36: = 


H- = x Ly Le 
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38 is identical to Eq. 5 withn H = F =0, which was from moment - 


material contained in > Appendix. 


and in the drain pipe it is 
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‘or the tank iquid and 
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OBLEM 


by Thomas J. 


low head radial gates writer is interested in the discussion 
_ The fourth method presented by ‘by Mr. Rhone appears to be identical to th the r 
procedure, used by the U. .«S. Army Engineer Waterways Experiment Station, to 
: develop a discharge coefficient chart for radial gates on curved spillway 
crests. — _ This chart was prepared for and included in “Hydraulic Design Cri- 
- teria” sponsored by the Office, Chf. of Engineers, Dept. of the Army. The 
WES chart (Fig. 1 herewith) and Fig. 6 of the author’s paper are quite simi- 
me except for omission of the data points on Fig.6. sal 
The bounda ry geometry of flow under a radial gate on a curved crest is 
complex. To simplify discharge computations many engineers have used gate 
openings and heads based on the gate seat and crest elevations. Furthermore, = 
in most instances, the effects of the direction of the approaching streamlines — 
_ have been ignored. Consequently, the results are often misleading when ap-— 
_ plied to other structures of f dissimilar geometry. 
“call The WES method involves the net gate opening between the gate lip and al 
curved crest. It also includes the head to the center of the gate opening as 
- suming theoretical orifice discharge and atmospheric pressure under the is- 
7 suing jet. Discharge coefficients were computed and plotted as a function of — 
_ the angle formed by the tangents to the gate-lip and to the curved crest at the 
4 net gate opening. It is believed that this angle controls the direction of the 
streamlines in the immediate vicinity of the gate opening. ‘The = should © 
therefore be applicable to geometrically similar structures. sa : ak 
_ The WES chart was developed using model and prototype po “The > data 
points shown are for spillway crests 3 shaped to the form of the lower nappe of 
a free jet over a sharp crested weir. _ The effects of gate opening- head ratios ‘9 
_ on the discharge coefficients were ‘studied. No well defined families of curves — 
could be drawn. se Therefore, an average curve was accepted for design ~ 
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DISCUSSION 
the fourth method for r the 3-ft gate opening is 4,880 cfs which is sent 11% 
a! greater than the measured prototype discharge of 4,400 cfs. The writer, in an 4 
attempt to discover the reason for this difference, computed a discharge of Fe 
4 ,482 cfs. The latter figure is within 2% of the measured discharge. A net ; 
gate 2.90 ft, an effective head of 30.75 ft, and a of 73.1° were de- 
termined graphically from available data. 3 A discharge coefficient of 0.681 _ 
from Fig . 1 was used in the computation. Thewriter would be interested in o 4 
seeing comparative figures used by the author for his computation. _ Sumer 
= for the 6-ft and 10- -ft gate openings were in close agreement 
a with discharges computed by Mr. . Rhone which are about 5% greater than the 
In the study of model it that the 
creased rapidly for gate openings less than 0.1 ft model. For similar proto- g 
type gate: openings the coefficients were essentially constant or showed 
slight decrease. This difference between model and prototype data is attrib-— 
uted to inherent difficulties in accurately measuring small gate openings and ‘ 
flows in models. Therefore, only model data with gate openings | greater than _ 
0. 1 ft were used in preparation of Fig. 1 
The writer concurs with Mr. Rhone’s comments concerning the need for 7 
- teolating the effects of such factors as approach depth, pier shape, gate radi- 
us, trunnion position, etc. However, it would also be necessary to determine 


i“ the effects of numerous combinations of these variables. A less elaborate ap- = 
be” ‘proach for early usable results would be to standardize some of the er, k 
TOCH,* M. ASCE. —As noted previously,> it that a 
general mathematical analysis ‘for the radial gate in a flat- -bottom, a 
_ lar section is not practical. Furthermore, it is felt that continued use of 
Horton’s article (3) asa ‘reference in the literature on Tainter gates is 
-grettable. Horton not only chose to equate a contraction coefficient toadis- 
Ra charge coefficient but did so for gate shapes of quite dissimilar geometries. me 
> Of the four methods presented by Mr. Rhone, only the last appears to on 2 


fer from the others. The first two formulas are evidently ‘the same, as ‘the a t 
constant value (2/3) 2g can be incorporated discharge coefficient 


cal factors in terms of the design 


Because of the algebraic difference first formulas and 


the fourth, the coefficient of discharge, for the latter, will have a quite dif- i . 
ferent numerical magnitude; it must still be en however, | on the same : 


: 


3 “Model- Prototype Conformance of Radial Gate Discharge Capacity, 
Dam, Montana,” U.S.B.R. , Hydr. Lab. Report No. Hyd 433, March 4, 1957. eh 
ae _ 4 Research Engr., Iowa Inst. of Hydr. Research, Iowa City, Iowa. oak 
= “Discharge Characteristics of a Tainter baa on a Spillway,’ a by 2. M. - Glowiak, 


8.1 Thesis, ciate Univ. of Iowa, 1955. = 
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ies: 
— 
= ore, to va rable. 
boundaries are first two only in the lip angle 
tics of the e iffers from that inclus >. 
The third form ionable to a nnarent that Fig — 
in@. It seems questiona! Also it 
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"variables c or parameters. In fact, the importance of this aspect of the problem _ 
can hardly be overemphasized: The coefficient of discharge must be expected = a 
tob be a function of the of gate and dam crest. In this 
5 ‘a connection with Mr. Rhone’s request for further information regarding» 
tests on radical gates, a study) made at the Iowa Institute of Hydraulic Re- 5 
’ _ search should be cited. Although this work did not include as many variables + 
= the author thinks necessary or desirable, Z. M. Glowiak did discuss and 
arrange the results he obtained in a rational manner. _ Another report® con- 
tains fairly complete data on flow conditions and pressure distributions ob- i. q 
tained on a model of a radial relatively high dam. 
— THOMAS J. ee, M. ASCE.—Messrs. Pariset and Michel have brought Fa 
q out some interesting facts in their discussion, and their —e of the co- 
efficient- -discharge curves” is, 


— as thin as the Edgard de Sonza Dam crest profile; in fact, the ip one 
free-flow coefficient was only 3.84. The larger values of the controlled flow 
coefficients, and the fact that for a given H the coefficient becomes essential-_ 

;: ly constant when @ is greater than 70°, rather than increasing, as shown on 
a Fig. 6, indicates that other variables should possibly be ‘considered in an ag 
= analysis. The use of the cavitation cocfficient presented is a very good meas- 
means of determining cavitation in conduits and entrances, ‘and in 
te the case of gate slots, has been proved with model prototype comparisons.9 . 
i The writer assumes that the statement made by Messrs. Pariset and Michel, 


t the irs. P = 


te is in error, and should read either “minimum allowable pressure,” or “ max- “fs 
mum allowable subatmospheric pressure,” as stated later, = 
Mr. . Bowman’ Ss awareness of the problems presented in the paper and his — 
"subsequent comprehensive | discussion are very gratifying. The increasing > 
 —— of papers that are appearing in technical publications seem to —— 


“cate a favorable trend in “sharing the knowledge gained from basic research. es 


It is hoped that this policy will continue, and will be as well received as in 
Bowman’s discussion of the various discharge quantity equations pro- 
ides a clear, concise analysis of the advantages and disadvantages of each. - 
_ He also provides a guide for the selection of the best combination of equations 
to use under various conditions of gate design, gate opening, and head a 
gate seal. His own method of obtaining an approximate rating curve is anex- 


Surmounted by Sector Gates), Calderini, Milan, Vol. 
- “Engineering Hydraulics,” edited by Hunter Rouse, John Wiley and Sons, p. 30. Sea 
“Hydraulic Characteristics of Gate Slots,” by J. W. Ball, Journal 
Proceedings, ASCE, Volume 10, October, 1953, pp. 


a | 
= 3 
: 
and con- — 
— 


-tractions, etc., a are, for the a of and will vary vary 

according to the practice of each design office. 

‘The datum profiles used in the series of t tests described in the paper were 
based on the hydraulic model investigations ; of the Hoover Dam spillway. 10 10 - 
_ Mr. Bowman’s discussion on gate seals is very informative. His observa 

: tions on prototype installations broaden the scope of the information provided | 3 

in the paper, which, as he stated, was principally derived from Bureau of 

Reclamation practices. _ The writer cannot answer the quasen as to how a 
L-seal negotiates the bottom corners of a radial gate. In most of the struc- 
tures investigated, the side seals were of the molded ciate type, and the bot- 
_ tom seals were usually of the rectangular type, with a tightly fitted butt joint | 


In his discussion, Mr. Cox has brought forth an unfortunate - omission; 
: one: that the writer did not identify the source of Fig. 6. _ The curves ' were . 
_ obtained from the Waterways Experiment Station chart, as he stated. - The 
_ discharge quantity that he questioned was given as 4 880 cfs in the table, but 
M4 should have been 4,580 cfs, which is within 4% of the measured —. y The 


writer’s figures in arriving -at this quantity aida as follows 


Net gate opening = 96 


the solution, and either could be accepted. Mr. Cox’s explanation for rome 
; larger values for the computed discharges at bigger gate openings i is entirely 
feasible. A steeper crest would permit higher flow quantities, other | conditions 
being equal, thana flatcrest. This is particularly true for the larger gate 
openings, as shown by Messrs. Pariset and Michel, and Bowman. The writer 
agrees with Mr. Cox’s statement concerning the standardization of some 
the geometrical design. Lc: The use of the design head as a basis for — 
dimensionless ratios is successfully used in other phases of spillway design, 
and should be equally satisfactory in radial gate design, neh 
Mr. -Toch implies that the discharge formulas presented in the paper poe " 
also recommended by the writer; this was not the writer’s spamegpead The = 
_ _ methods presented in the paper are inmore or less common use, ‘and the pur- 
‘ - pose in analyzing and discussing them was to show that each method is inade- 
quate to some degree. The similarity between the first two formulas is in the 
_ form of the equation only. It should be obvious that the coefficient-of-dis- a 
4 charge values in the first equation are average values obtained from numer- = 
ous model studies, irrespective of gate geometry. In the second equation, the — } 
coefficients are related to the geometrical configuration of the gate. Appar- oa 
eatiy, the originator of the third method believed that the use of any gate ona 
crest would affect the coefficient of discharge to the extent shown in the lower 
graph on Fig. and that the gate geometry would further ‘modify the coeffi- 
_ cient in the manner shown in the upper graph . The writer cannot agree that , a 
the upper: merely relates the cosecant and cosine of any given angle,” 


“studies Crests. for Overfall Dams, Boulder Projects, Final Reports, 
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stated by Mr. a 
simplified method for obtaining the value of 1/sin by using the readily 
available values of the gate radius, trunnion height, and gate opening above 4 
crest. Mr. Toch— is correct in believing that the coefficient of discharge is 
also. a function of the crest shape. . Previous discussions have vividly demon-— ; 
strated this to be true, and it should an important consideration 


| 

| reviewed 

% d _ The thesis by Mr. work 
Toch ar d it is a comp discharge charact 


_ DISCUSSION OF RELATIONSHIPS BETWEEN 


Closure by Walter L. Moore 


| WALTER L. MOORE,! 


th the inconvenience of. “the Colebrook- -White relations for resistance 
problems, presented different approaches to aid in the achiev ement of. a ‘more a 
_ convenient method. — Emphasis was given to the limitations of the Colebrook- 7 


will be cleared up to the extent 
Mr. Ackers’ discussion was particularly interesting in that he described a 
different attempt to present the Colebrook-White relations ina more ‘usable 
form than previously available. The first of the two papers to which he refers — 
“gives an excellent comprehensive review of the pipe-resistance problem. 
essence, 2, thea approach described by Mr. Ackers consists of plotting the Cole- 
-brook- White relations in terms of dimensionless velocity, size, and slope 
& parameters. It is of interest tosee how the parameters he used can be formed \ 
by the combination o of the more conventional dimensionless parameters of 
friction factor f, Reynolds" Number Nr, and relative roughness e/D. Ackers 
dimensionless parameters are listed below, together with their definition from 


hi h l 
is paper and elr equival ent in terms | of Conventional 


-1n2 (5 ) 
Ackers’ Eq. 1 the k is equivalent to e as used in the rest ‘of the 


om _ The general resistance diagram, plotted in terms of the foregoing para- 
_ meters, represents a logical arrangement and is based on sound dimensional 


involving a simple single pipe. pram direct solution may be obtained for either . 
the velocity, the diameter, or the head loss » if the other necessary quantities 
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Ss. . It has ‘certain advantages which are readily apparent for problems 


— 
amore comprehensive presentation of the resistance problem. Some of the 
_ 
. 
— = 


1980 
problems i in which the diameter, the velocity, or the discharge are unknown. _ 
_ ‘The new diagram proposed by Mr. Ackers is not an improvement in all re- 
_ spects however. Since each of the parameters V, R, and S varies over a wide 
range of practical values all of the scales must be rather compressed, tn 
_ impairing the precision with which values can be read from the charts. ot 
course, large charts can be made to. improve the precision (see charts in ol 
pouch of Ackers’ reference 7) but aan become rather awkward to use and 
are not to be widely available. On the her hand, the pipe- -resistance 


Moody diagram. the ‘scale for f can be chosen chosen large enough to read 
The physical : significance of the variables { is 5 not as readily apparent in the ; 
 Ackers’ diagram as it is in the Moody diagram. .. The interrelationship c of the 
velocity, size, and head loss (or slope) is apparent for the proposed diagram, od 7 
7 but don is already apparent from the general form of the various resistance 
The significance of changes in the ‘Toughness ¢ or the fluid viscosity 
are pial readily apparent from the new diagram, because these variables oc-_ 
cur in more than one of the dimensionless parameters. _ The Moody diagram, 
on ‘the other hand, clarifies the significance of these variables and brings o bal 
- more clearly the meaning of the “smooth pipe” and “ rough | pipe” limits. i. 
Fer complex pipe systems, “including transition losses or pipes in series, — 
ae _ the writer can see little advantage in the Ackers diagram as compared with 
the Moody diagram. If the rate of discharge or the velocity is known, either 
a diagram may be easily used. If the discharge is not known a trial solution will - 
be required with either of the diagrams. = 
ow? writer has found the exponential relations so useful. The writer believes itis 
- more rational to examine the range of Reynolds Number and relative rough- | 
e ness likely to be encountered and to develop a suitable exponential resistance - 
- equation ‘from the modified Moody diagram, according to the methods he has 
i _ proposed, than to arbitrarily “select” an exponential formula without due con- 
- sideration for the range of parameters involved in the particular problem. 
Mr. Szesztay briefly reviews applying the Colebrook- White 
formula to a graphic ‘solution of channel flow. Essentially his method consists 
of specifying the shape parameters of the cross section in terms of the water 
depth: h, thus ‘making the area and hydraulic radius a function of water depthh | 
only. Fora a given sand- roughness k the discharge may then be ex- 
' pressed as a function of water depth h and ‘slope by means of the Crump 
_ the Colebrook- White formila, as given in Szesztay’ s Eqs. 5 Ay 
and 4, The chart presented in Fig. 3 permits a ready solution of these rela-— ¥ 
tions. One could easilyalter the relations given by Mr. Szesztay for the chan- 
; nel shape as a function of water depth, if he wished to use different criteria : 
7 _ Although Fig. 3 .3 is clear for for the assumed k-value, it is not clear to the 
; writer how the corresponding va value of n = 0.025 can be ‘stated. tt would seem | 
q that the value of n corresponding to agiven roughness value k would bea func- — 
tion of the hydraulic radius and, hence, of the water depth a: Seine : 
_ States that statistical studies of n and k apparently fail to reveal any con- 
sistent effect of the value of the hydraulic radius R. It would seem that this 7 


It is for these complex pipe systems, well as el pipe | networks, that 
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In Fig. 3, the key for the sl slope S is. listed as witiaiinee over centimeters ..' 
whereas it appears that it should be ‘centimeters over kilometers. 
_ The writer was pleased to note the similarity of viewpoint presented by ; 
Mr. Miles. His extension of the concepts to a tabular comparison of e/D and 
Ha; en and Williams’ C-values was particularly valuable. The results that the 
-Hazen-Williams C approaches 160 as the relative roughness approaches 0 was 
interesting. The reasonableness of this value is indicated by the C of 155 for — 
new bituminous enameled- lined pipe of 16 in. or larger pee in the 5th edition 
Babbitt and Doland’s “Water Supply Engineering. 
a! The writer is grateful to Mr. Bilonok and Mr. Soiiiininiins for correcting a — 
numerical error in the first example of the appendix. _ The figures shown give 
: a value for "hg of 4.88 ft rather than 5.84 ft, as printed in the paper. . The cor- 
rect figure results in a value of Ky = 0. 455 thus giving a nn 
for this s example of: 
: This aeeniiiiaa erroralso led to an incorrect value for the lengthof the 12-in. 


Z _ steel pipe appearing in Table 2. In order that the , equation for the 12-in. _— 
pipe, hey = 3 1.99 Q1.88 be correct, ‘the length of the pipe should be be 4,370 ft in 


place of the 3,660 ft shown in Table 


Mr Bilonok also recognized that are in common use 
by practicing engineers because of their convenience, as compared with the — 
general resistance diagram and the Darcy-Weisbach formula, : 
In discussing the effect of aging on pipe resistance, Mr. ‘Bilonok states that — 

_ roughness factor e, is a “complete measure of the size, shape, and distribu- 
tion of individual manees elements.” Although the introduction of roughness ? 

_ parameter e indicates a great improvement over other methods of charac- — ; 
_terizing surface roughness, the writer agrees with Mr. McPherson and Mr. 
Ackers, that shape and distribution of individual roughness elements are not 
“necessarily characterized by the value. of e. For the case of statistically as- 7 


treme variations of such as grooved and corrugated 
some other geometric parameters” are important | besides -& measure of the | 
‘The writer is puzzled by Mr. Bilonok’s question ‘regarding the finding of 
exponent in 2 where m = 1.95 » which he says, “ is not in agree-_ 


Nr and f, c or r e/D listed ; are e used to enter Fig. de of the paper, the value | of m 
will be found to lie between 1.95and1.96. 
.: The writer is indebted to Mr. McPherson for a very painstaking and com- 
_ prehensive discussion which stressed some points which probably deserved _ ; 
4 more emphasis. Also, Mr. McPherson copnyertly: misunderstood a number of 
io It is obviously true that no system for approximating a complex relation- Da 
ship can be more precise than the relationship being approximated. _ Thus the 
proposed exponential equations are certainly no more precise (but fortunately — 
_ only slightly less precise within the limitations imposed) than the Colebrook- 2) r 
White e relationship, from which they were derived. Mr, McPherson states that 
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values ¢ extend beyond the range of water- -works "practice. It was not 
7 writer’s intention that the concept be limited only to water-works practice. ce 
SS Mr. ‘McPherson questions th the statement that the exponential equation will 
: i agree with the ger general resistance diagram \ within 5% for a 20-fold variation S| 
Q. This statement is significant, not because one expects to analyze a water — a 
_ distribution system over a 20-fold variation in Q, but in that it indicates that 7" 
a . satisfactory exponential equation may be developed « even when the initially 
"expected value of Q is poorly chosen. It must, of course, be realized that for 
the exponential approximation to apply | for a full 20-fold variation, it must be 
developed about a central point in the range rather than about one extreme > 
end of the range. Thus, in the first example in the appendix, which gives a 
head loss of 4.88 ft at a Q of 3.53 cfs, , the comparison should be made at | 


os at Q = 3.53 V20 = 15.8 er 


From the equation} 0.455 88 


= 0. 290 ft at = = 0. 79 
ft atQ = 


percentage di differences are as follows: 
= 0.79 cfs; (0. 292/0. 299 . - = 2% 
The writer agrees with Mr. McPherson that in neetiatian of water-works | net-— 
works a constant value of the exponent m may normally be used. The writer : 
is more rational and will also | give better “results than arbitrarily using 
1.85 exponent from the Hazen- Williams formula. Mr. McPherson’s discussion 
: does not seem to recognize that network problems also exist in fields other — 
than water-works practice. In some instances it may even be desirable to use 
different exponents m for different branches of the system. Contrary to Mr. 


McPherson’s digital computers are not necessarily restricted to 


“the preciseness of the general resistance diagram implied the ‘ exponen- 

tial equations’ of the author does not appear to be justified.” It is a misunder- 
standing of the original paper to consider that the exponential equations imply 
about the preciseness of the general resistance diagram. They are merely an 
— attempt to put the information in the general resistance diagram in a more | . 3 Re: 
convenient form for a number of types of resistance problems. .x+_ 
McPherson is concerned because the numerical values in Table 2 pass 
the realm of water-works design practicality, 4 
[a jj The table was designed merely as a means of illustrating the effect of a [i 
j = 
4 a 
From the general resistance diagram a 

= 0,299 tat Q = 0.79 chs 0.019) 
“Ao, 
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= exponential resistance equations,” which was defined in the 
: mediately following Eq. 2. Thus, any equation of the form hg = kK Qu (which 
ae the Hazen- Williams formula asa special case) is considered an ex- 
ponential formula. — Mr. McPherson must have misunderstood the term ine 
4 order to make the statement that he was “forced to reject as. - obviously un- 
realistic (and clearly not proven) the statement by the author: ‘although the 
methods for analysis of networks will not be considered here, ‘it is pertinent _ 
to note that exponential resistance relations will be most convenient for use in - 
= type of analysis.’ McPherson goes on to propose the use of the 
Hazen-Williams for network analysis which is really a confirmation 
of the writer’ s statement although restricting i it to the ‘special case of an ex- 


Again Mr. aiiantiei apparently misunderstood the statement that, “thus 
ie Hazen-Williams formula with its exponent of 1.85 corresponds to a com- 
a clearly deals with the effect of m on the change of head loss with 
change in Q. ‘The statement should thus be taken to indicate that only along the 
m = 1. 85 will the Hazen- -Williams "equation have the correct exponent. 
_ Obviously the whole field of friction factors and Reynolds numbers may be ~ 
_ covered by lines corresponding to an exponent m = 1. 85. _ This is what is done 
in the chart presented by Mr. McPherson as Fig. D- 2. The important point is A 
ze the change in f with Reynolds number occurs at the same rate for both 
the Hazen-Williams equation and the general resistance diagram, only along ~ 
the line corresponding to m = 1.85. This may be clearly shown by plotting» 
coordinates of the line m= = 1.85 from the modified Moody diagram onto Fig. 


representing constant relative roughness on the general resistance diagram. 
To | the ‘right of the chart the Hazen-Williams equation would indicate a con- 
tinuing decrease of friction factor with increasing Reynolds number, whereas 
& rough pipe relations represented on the resistance diagram show that the 
friction factor should become constant. It might be argued that C. H.w, could 5 
: be decreased as Reynolds number is increased, thus forcing Seana Williams 7 
equation to yield a constant f at high Reynolds number. This would complete- 
ly destroy the added convenience of the a _type equation in this" 

range. 

‘exponent m = 2 in in order to fit the conditions’ for a constant f-value. rte 
Further evidence of the correctness of the statement is given by 


Hazen- Formula. ” This equation 


because D0-015 will be. 


- practically 1 for any reasonable value of D. Thus if Cy is set equal to 
— 1060/Cy w, and and m = 1.85. , Eq. 5 will be se seen to correspond to the foregoing — 


“Relaxation Methods for Pipe Networks,” by T. Chapman. Civ. and Pub. 
‘Works Review, Vol. 51, No. 605, 
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data by in “Tables D- D-3 and D- D-4 for new 

_ linings and old linings is of interest. If the C-values given for Reynolds num- 

7 _ ber and friction factor are plotted on the modified diagram, the points for the | 
new linings, for the most part, fall between the lines of m = 1.90 and m = 1.96. 
For the old linings (Table D- 4) the points, for the most part, lie between the - 

lines m = 1.94 and m m = 1. 975. This data, also, indicates that the resistance — 
relations for the older pipes would be better. approximated by larger 

: of the exponent m than for the new pipes. It is also of interest that the indi 
cated values of m are all considerably about the 1.85 line corresponding ~ 
= Hazen- Williams equation. a Thus this data also supports the writer's con- 


approaching 2, would be preferable for old water mains. 
on It is the writer's belief that the modified Moody diagram in conjunction 
with available data for e, or better still test results such as presented by Mr. 
- McPherson, will enable the engineer to select the exponential equation which * 
_ will be most appropriate for the type of problem with which he is dealing. — d 
* a writer recognizes that there is a possibility (which at this time seems _ 
7 ¢ rather remote) that large changes in roughness form due to encrustation and 
the forming of ridges or waves might alter the form of the transition function r 
from “ smooth” to “ rough” pipe conditions and thus affect the proper value of a 
‘the exponent m. That the shape of the ‘roughness : should change enough to —_ 
__ nificantly affect the treatment function seems unlikely, and in the absence of © 
_ reliable experimental data on this point, the modified Moody we appears 
to be a reasonable improvement over existing methods. by 
Additional research to explore the shape of the transition function and the 
- geometric form of roughness for encrusted pipes would certainly be of inter- 
est. In such an investigation it may well be that the reduction in diameter will 


into the ‘equations to the fifth power 


The work described by Dawsey attempting to relate measurements of oii. 
ness profiles to the equivalent sand grain roughness e has much appeal. a ; 
be able to run a test on a small portion of a conduit lining, or better still on a 
- cast taken of a portion of the lining, and convert this to a hydraulic roughness 
oe parameter would be a tremendous forward step. _ There are many difficult 
he problems of measurement, and meaningful analysis of the measurements, that 
- will have to be solved before this approach is likely to yield reliable results. 
ne, The magnitude of the potential gains are so great that efforts toward this goal 
seem well directed. The writer thanks Mr. Dawsey for “noting | the omission of 
parentheses in ‘the | ‘Colebrook-White “equation as presented in Table 1 of the 


The writer is grateful to all of the discussors for contributii 4 ae 
__ information and new thoughts regarding pipe-resistance relations. The dis- 


cussions have supplemented the paper by presenting other attempts to _ 7 
velop convenient as well as reliable methods of evaluating pipe-resistance. 
: is hoped that the concepts presented will lead to an improved understanding « 2 


a pipe resistance ‘relations. Time and the experience of the practicing enginee 

_ will determine which of the various — will prove to be most useful, 
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CHANNEL-SL SLOPE FACTOR IN FLOOD- FREQUENCY ANALYSIS 


MANUEL A. BENSON, M. ASCE.—Mr. Sammons is to be commended for 


his his thorough analysis and reanalysis of the writer’s data. The original article — 
7 was intended to cover only the choice of a simple, yet efficient, index of main-— 
channel slope. Mr. Sammons covered a wide scope, including basic data col- 
lection, mapping standards, flood- -frequency analysis, multiple en, 
statistics, and the use of automatic computers. Much of his discussioncon- 


. cerns points which will be explained in detail in the pending final report of the © 


{,, ‘completed project. ¢ Mr. Sammons makes some sweeping indictments of vari- 
ous practices, to most of which the writer wishes to plead “not guilty.” _ For 
example, the writer sincerely believes he does not share that misconception 
q of the use of automatic computers, and he is certain ‘that he does not apply 
Mr. Sammons puts much stress on the distribution of the variables which 
used in the regression analysis. Multiple- -correlation techniques, to be 
valid, require linear relationships between the variables. ‘The transformation 
to logarithms of the original variables was made on the basis of comparative | 


logarithmic and plotting between the dependent and independent 


independent variables. variable (surface- storage area in ponds) ‘used in 

the final analysis was adjusted by a “ rectification constant” determined graph- 
__ ically to be necessary to linearize the logarithms of the variable. yen (as > 
= by Mr. Sammons) a 

eventually used were plotted on log-probability paper and were proved to be 

 emuitialls linear, within 2 standard deviations of their means. This. con- 
the validity of the logarithmic transformations. 

No logical reason is apparent for Mr. Sammons’ statement that the curvi-— 
linear: regression constant of “a” in Fig. 3, indicate that “a = 


_ logarithmic model would best fit the data.’ ” Elsewhere, Mr. Sammons refers — 
to the “ concave upward bias | of the constant term a on Fig. 3.” A bias is much 
os apt to be introduced by an a priori theory than by adherence to the data. 
The distribution of most of the variables, which may influence flood peaks, 
is ‘sufficiently close to log-normal, thus that distribution may be assumed. In 
any event, _ their exact distribution is not critical for the purpose of relating 
- them to flood peaks. _ However, the writer is not willing to assume log-normal | qj 
distribution for the annual peak discharges at a station, because this is a 


cxttical assumptio n. The form of the distribution becomes the principal basis q 


7a 1 1 Hydr. Engr. Research, U.S. Geol. Survey, Washington, D.C. 7 
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for determining extreme flood peaks. Statistical me methods may be nonpara- 
¥ metric as well as parametric. _ The writer prefers to _correlate topographic 
climatic variables with the “percentiles” of the flood-peak distribution 
_ rather than with its parameters (mean, standard deviation, and skew), simply a 
because the use of parameters means making assumptions: that not > 
necessarily or generally true and that may actually be in error. 
P Mr. Sammons cannot suppose that the investigation which he cites, i in 
which precipitation proved to be more important than the drainage area, is” 
= applicable or reflects on the results for New England. If the analy- 
sis were made for a group of small drainage areas (necessarily limited in — 
: } range of sizes) within a region of widely varying precipitation, it is to be ex-— 
pected that precipitation will: prove more important. In fact, if all the drain- 
age ‘areas are the same size, drainage area would be of no importance what- 7 
soever, at least as shown by the multiple- regression analysis. = a? 
_ The writer agrees with Mr. Sammons that the median values of variables 
: might be better indices than the means. _ However, the ‘U.S. Weather ‘Bureau 
publishes only the means for climatic factors such a as precipitation and tem- 


perature. _ The mean of topographic variables, such as elevation and land 
slope, is more simple to compute than the median. ee 

scales of topographic maps New England were adequate to define 

- profile elevations for the basins used in the analysis. _ Regions of flat slopes, 
ef less adequately mapped regions, or extremely small watersheds might re- 
: quire surveys or photogrammetry to determine the channel slope. The writer 
did not investigate the effect of variation in map scales because the available 
introduced no problems forthe data studied. 


_ The related flood peaks of recurrence intervals, 


2. main-channel slope; 
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ROLL WAVES AND FLOWS IN OPEN CHANNELS? 


Ishihara, Yuichi hi Twagakt and Yoshiaki 


TOJIRO ISHIHARA, 1 YUICHI IWAGAKI,2 M. and YOSHIAKI 
IWASA. 3_The writers have reviewed this paper with interest, because we 
also have made continuous investigations of the hydraulic characteristics of z a 
roll waves 2 and other associated problems in laminar and turbulent flows from 


= fortunately, peo our results are quite different from the author’s ‘study, 
= 4 and we also. do not follow his treatment of roll waves and slug flows from the 


knowledge we obtained through research in the hydraulics laboratory of Kyoto 


ss waves, needed by the engineering request for soil conservation, have been a 
_ studied by the writers 7 others since 1950. Several reports were already ’ 
— published in Japanese*-7 and in English8, 9 after completing specific problems 
to roll waves. ‘Through our research works, the basic concept of the 
- treatment for the hydraulic characteristics of roll waves is essentially due to 
the studies of H. Thomas!10 and R.F. Dressler.11 More precisely, the deriva-_ 
7 tion of the hydraulic characteristics of roll waves is based on the fact that 


1 Dean, Faculty of Engrg., and Prof. of Hydr. Engrg., Dept. of Civ. Engrg., Faculty 
Asst. Prof. of Hydr., Prevention Research Inst., Kyoto Univ., Kyoto, 
Asst. Prof. of ‘Hydr., of Engrg., Faculty of of Engrg., , Kyoto Univ., Kyoto, 
> 
ane 4 “Theory: of the Roll-Wave Pee Laminar Water Flow on a Steep Slope Surface,” i 


by T. Ishihara, Y. Iwagaki, and Y. Iwasa. Studies on the Thin Sheet Flow, 5th Report.. oa 
Trans. JSCE, No. 19, April, 1954, (in Japanese), 
5 “On the Hydraulic Characteristics of the Roll- Wave Trains,” by Y. Iwagaki and 
wb Studies on the Thin Sheet Flow, | 7th Report, Journal JSCE, Vol. 40, No. 1, 
ee “Mechanism of Water Erosion on Land-Surfaces and Roll-Wave Trains,” by T. 
Ishihara, Y. Iwagaki and Y. Iwasa. Bulletin, Engrg. Research Inst., Kyoto Univ., Vol. 7, | 
March, 1955, (in Japanese). dt ain Wao Flo? by 
" “Fundamental Studies on Land Erosions due to Rain Water Flow,” by Y. Iwagaki. 
Dr. Engrg. Thesis, Kyoto Univ., March, 1956, (in Japanese). = 
ak “On the Rollwave-Trains Appearing in the Water Flow ona_ Steep Slope Surface,” _ 
_ iby T. Ishihara, Y. Iwagaki and Y. Ishihara. Memoirs, Faculty of Engrg., Kyoto Univ., 
2 9 “The Criterion for Instability of Steady Uniform Flows in Open Channels,” 
‘Iwasa. ‘Memoirs, Faculty of Engrg., Kyoto Univ., Vol. 16, No. 4, October, 1954. 
_ 10 “The Propagation of Waves in Steep Prismatic Conduits,” by H. Thomas. Pies! 
‘Hydr. Conf., Univ. of Iowa Studies in Engrg., Bulletin 20, 1940 
cs ri 11 “Mathematical Solution of the Problem of Roll-Waves in Inclined Open Channels " an 
by R.F. Dressler. Comm. on Pure and Vol. 2, No. 2/3, 1949. 
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(1) the wave velocity of roll waves is definitely constant (terminal velocity by 
2 author), regardless of their flow regime, when the final shape in wave — 

_ pattern is obtained, and (2) the flow velocity is, in all cases, less than the 7 
wave velocity of roll waves in magnitude, thus the flow regime | of rear side of 
a roll wave is shooting, whereas that of the front ‘side, including the wave 
crest, is tranquil when the train of roll waves is observed from the moving 
coordinate system travelling at aconstant speed of wave velocity. The former 
condition was also verified by the author, , while the latter condition of Thomas 
was not equivalent to the author’ s in all points aaremnenee? the author divided he 


_In order to ascertain our standpoint cilities the author’s a for 
7 the hydraulic characteristics of roll waves and slug flows, our analytical 
Z rs treatment o on the hydraulic characteristics of roll waves and t the experimental : 
4 verification by the direct measurement of roll waves by means of electronic aa 
a _ Based on the principle of a constant wave velocity indicated in the fore- 
7 going, | the unsteady behaviors of roll waves can be reduced to those in steady 
ye regime by eliminating the time derivatives in the basic quasi- linear partical | 
_ differential equations of hyperbolic type for unsteady flows, when viewing the 
| wave pattern f from the moving origin. 7 The hydraulic characteristics of ‘roll ; 
waves are then easily derived in the following. 
a The equation of momentum for unsteady flows in apes ‘chante with | a con- 


= 
in which the pressure distribution is assumed hydrostatic, x is the distance — 
from the « origin n along the channel bed, t represents time, is the mean 
Velocity of flow, D denotes the depth of water, A is the flow area, R nhs ronal 
_the hydraulic radius, C ‘is Chezy’ s coefficient (assumed variable), g denotes — 
the: acceleration of gravity, @ is the inclination angle of channel bed, anda — 
oS “equals the momentum correction factor usually defined. Apparently, the _ 
= profile of laminar flow is of the parabolic type which was already de- . 4 
4 rived | by many hydraulic engineers as well as the author and the momentum 
correction factor becomes then a constant of 1.2 with a simple calculation. 
the other hand, the evaluation of a in turbulent flow cannot be made in a defi- 


nite number, because it is still a function of flow and channel 


A but it is satisfactorily assumed 1.05 in engineering practice. In the derivation 
= of Eq. 1, the curvature of a stream line is neglected, and thus the influence of — 
5 ‘surface tensions is also ignored, which is different from the author’s treat-_ 


Taking account of the preceding principle for a constant wave velocity, the 
_time derivatives in Eqs. 1 and 2 can be ttminated by setting { =x - Ve t, 
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: In reality, the wave velocity may not be considered constant | in a channel of 
irregular section, but the treatment of analysis will be proceeded from a ~ 
macroscopic aspect for mean flows. Eliminating V from Eqs. 3 and 4, the 2 
_ steadied wave min observed from the moving origin can be renee 


in which K isa progressive ¢ iasdneai rate defined by (Vy - V)A=a constant. 

_ As seen in Eq. 5, the wave pattern of a ‘roll wave cannot be derived in a peri- a 
_ odic form observed in laboratories, and thus the solution for roll waves must 

be obtained as a discontinuous wave pattern by combining each of water sur- 


profiles the shock condition as Dressler did . Dressler a also ob- 


The actual wave pattern, however, observed in channels is rather in a similar 


to discontinuous periodic s solutions. ‘ 


file, as the flow regime changes from tranquil to shooting, so that in Eq. cdl 
both numerator and denominator become simultaneously zero and all the hy- 
draulic characteristics of roll waves can be uniquely determined at this sec- = 
tion. The resulting expressions for the mean velocity and the ey radius it 


= 


a(dA/dD) - Vala- 1)(aA/dD)2 +{ (dA7aD) 2. 


a(dA/dD), - (So/N2, 
in is the » jetted perimeter, Ne denotes the Froude number defined b by 
Pg Vv o/ (g Ro cos g)i/ , and the subscript o indicates the values at the control | 
- section which are eventually equal to the quantities in uniform flow. The p pri 3 
. A gressive discharge rate can also be calculated by the use of Eqs. 6 and 7. b Mealy 


the flow is laminar 1 in a channel, the normal mean 
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indicated as of ‘the wave veloc ity and Froude number as 


2 


7 Water-surface ‘profiles as each wave pattern of roll waves are obtained oo : 


integrating Eq. 5 upstream and downstream from the control section. Al 
though water- -surface-profile equations cannot be integrated analytically in 
, those for laminar flows and for turbulent flows, characterized a 
by a constant. Chézy coefficient in rectangular channels, are integrable. In | 
cases of two-dimensional laminar flows, the resulting water- profile 
_ is indicated in an an inverse functional nal form as follows: coe 
A 


q 


in the constant is selected so that D = Do a - = and 
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Discussion 


is expressed by the three conditions: 


om rate ont ten . 14 and 18, th the e relationship between the the lowest and and high- nm 
est depths of a wave becomes 


(17 )} = cos (Dg - Dgp) - 1) 
fi = Cos - p(Dp - Dgp)} 17) 


Denoting the wave length of a wave L, the following relat relationship must must 


exist between two successive waves: 
— the water-surface profile for a particular roll wave is obtained by inte- 
ee numerically Eq. as complete solution of discontinuous gad 


| DS cos 6+ + (6/5)(K2Dj/e D2 Do) 


loge 


cos + (6/5)( (KDp/g D DB) + + (6/5)(K2/g Do) Do) DB 


De) 


j. The complete discontinuous solution Of To ves obtained bycom- 
pape ; _bining each of water-surface profiles through the shock condition, as will be (i 
— 
i | 
— 
a the subscripts of f and b indICate te values at tie lOWest alu 


‘The criterion for roll waves to ‘maintain their final wave pattern, 

_ seribed in the foregoing, v will be studied next. It is of common > observation > 

- that the wave profile of a roll wave is monoclinally 1 increasing from the rear _ . 
- side to the front side, so that the aolheting relationship for the water depth ue 
be obtained the when the waves are observed from 


With t the aid of Eq. 5, the above condition poe mene ered characterized by ; 


or, expressing in terms of Froude number, 


= 
increases to unity with increase in width of a channel. 
If the flow ‘turbulent, — = calculated in the same man 


ner, and it is 


for the channels i in n which the flow area is proportional to a power of the water Bs 


1+b ‘sin 


acteristics, 80 > that the Chezy’s coefficient C 
Inserting Eq. 26 Eq. , the condition | becomes 
+ b)mVo) 


i is s evidently the Vedernikoy nun number for the initial instability of shoot 


the velocity formula of Vedernikov’s type will be introduced as 

— 

¥, 

4 

| 

| if 
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and | in 1 specified c ‘cases of particular ve velocity formulas 
Jeffreys14 and and Patterson19 obtained the instability conditions. 


roll waves frome in a their final wave pattern, is ; that two conditions | for main- 


x 


The previous e: explanation oor the hydraulic characteristics of roll waves is i 3 
a a summary of our past analytical treatment developed by means of the basic ~— 
concepts of Thomas and Dressler. The experimental works for the charac- aT 1 
teristics of roll waves, to clarify the hydraulic behaviors of waves in definite 4 
forms and to verify the theoretical deduction: wave characteristics, 


gE — ject to laminar roll waves , the aluminum channel of. 20 ft in length was used, 
i whereas the experimental runs for turbulent roll waves, which are called as 
slug flows by the author, were made in the lucite flume of 36 ft in length and | a 
& of variable slopes from horizontal to 30°. Continuous recordings of wave pro- 
a files of roll waves and of wave velocities were made by the electronic. wave * 
_ recorders of resistance type, and the discharge measurement was made by _ 


volume metering of water at the downstream end. The details of experimental 


_ procedures were reported in the writers’ paper.* F Comparing the theoretical _ 
results described in the foregoing, the experimental behaviors of roll waves oe 


obtained by the writers will be presented in the following. | Data = he author > 
and Binnie16 will: be also supplemented. 


_ For the sake of simplicity in illustration, ‘the flow is assumed Sal 


"treated first. 
two- dimensional. In this case, ‘Eq. 23, the for lami- 


Re “The Criterion for the Possibility of Roll Wave Formation, A. . Craya. Proc. 

13 “Resistance Effects on Hydraulic Instability,” by R.F. Dressler and F.V. Pohle. | 

Comm. on Pure and Applied Math., Vol. 6, No. 1, 1953. 
_ 14 “The Flowof Water in an Inclined Channel of Rectangular Sectio on, ” by H. Jeffreys. 7 
Phil. Mag.,Series 6, Vol. 49,1925. = 

15 “A Criterion for Instability of Flow in Steep Channels,” by G. H. Keulegan ce 
G. W. Patterson. Trans. AGU, Part 11,1940. 
“Instability in a Slightly Inclined Water .M. Binnie. 


Fluid Mech., Vol. 5, Part 4, , May, 1959. ; 


the same expression of Vedernikov, regardless of essentially different ap- 
deduction. The formation of roll waves, therefore, will result 
‘thee instability of open-channel flows as many hydraulic engineers de- 
University. test tlumes Of a rectangular section were used ior this pur-— 
nose and in the earlier etage of rec eh nrocrescs which was mainly suh- 
— 
| 
her given bv V. D./ v _and the same rri- 
4 
4 


wt consequently, it depends on the velocity formulas like Chézy’ sor Man- 
-ning’s and the momentum correction factor. Since the Froude number is ex- 
pressed as a function of the Reynolds number, Manning’s roughness coeffi- 

_ Cient, slope and kinematic viscosity, when the Manning formula is used, the a 
instability criterion can b be graphically illustrated in terms of the slope and a 
Reynolds number as parametric expressions of the roughness and the mo- — 
mentum correction factor. Fig. 1 illustrates the criterion for instability of = 


flows and formation of roll waves obtained in the preceding manner under the 


_ assumption that the kinematic viscosity is 0.01 sq cm per sec. Two zones a 


Stable and unstable flows will be apparently observed in the figure. In open- | 
channel flow, the transition from laminar to turbulent flows takes place a 
Reynolds number from 500 00 to 1,500, - as often nm indicated by many engi 


Turbulent 
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Non-Formation | & Roll Waves , Mayer ~ Well Developed, writers 


Roll Waves Confused Pattern, @ Developed, 


The ee of the writers! found that the itis formula was not suitable in 
the transition region, but (1 + + b)/a was commonly very small and it became 


_ negative in particular cases depending on the channel slope, a the ex- 


oe Now, let us consider the behavior of roll | waves in formation when the flow 

oo is carried in an open channel. a For clarification in n explanation, the channel 
slope is assumed 0.02. -When the regulating valve for discharge is opened 
| from closure, the flow begins to be carried. The roll waves are not observed — 
Until the Reynolds r number becomes 50 50, as previously i indicated. With in- 

crease in discharge, the laminar roll w waves become appreciable. the 
‘Reynolds number increases to 500, so the transition from laminar to turbu- 
lent flows begins to take placeand the roll waves disappear. This phenomenon 

i -may be corresponding to that called the confused pattern by the author. | With 


the further increase in the flow tarbulent, and 
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again the roll waves are are predominant in forms of bore if the eo & j 


“ness is approximately less than 0.01 in Manning’s roughness coefficient, 
whereas for rough channels in which the coefficient is larger than 0. 01, the 
waves can not be observed. When the discharge is still increased, the roll 
; waves ultimately appear even in such rough channels. This fact is evidently 
- equal to the results obtained by Thomas, 10 Dressler!1 and others. = = 
4 ins In Fig. 1, our experimental data for roll waves are plotted, adding suppl 


4 


mental data of the author and Binnie. It is then concluded that the foregoing 
_ treatment for the formation of roll waves is verified by the experimental ob- 


servations. 2 Evidently, the law of resistance and the momentum correction 
_ factor play a very important role in n the criterion for formation of roll waves. 
Fig. 2 indicates the criterion ‘in terms of saemasiaed formula and momentum ~ 
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‘FIG. 2.- —ILLUSTRATION OF CRITERION FOR FORMATION OF ROLL WAVES It IN 
- TERMS OF LAW OF RESISTANCE AND MOMENTUM CORRECTION FACTOR. — 


the transit factor, and it _ illustrates that roll waves “= be observed in 


channels will be obtained as follows. 


in which is a a half of the angle two side 3 the 


same criterion for triangular channels, and it is obvious that the roll- -wave ? 
_ formation in channels be compared with that in rectangular 


—: behavior of wave velocity of ey waves when the the wave sation becomes 


ulti | mately will the use | of ons for 


4 
— 
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=Vy rectangular channels, in which do is 


As the dimensionless progressive discharge rate RY: in wh 


= K = 


- Vala - 1) + 24 
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is a function of and Np. ‘so the wave also ex- 


pressed | in terme of Froude number. 
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—CRITERION FOR FORMATION OF ROLL WAVES LAR 
é Fig. 4 shows the behaviors of Vw" for both laminar and turbulent flows. In the 
same figure, our experimental data are plotted with data of the author and 
Binnie, and it is found that the verification, to the theoretical 


the similar procedures, the behaviors of wave length, period, and 


as the 1e hydraulic characteristics are the 


a. progressive discharge rate and for the mean velocity of flow at the control 4 y s 
section, the wav ae 
— — 
q 
7 hich 
- 
(1/N#) 
: 
7 
— 
7& 


ollowing d dimensionless parameters of wave length, - wave parted and wave 


(g cos ), 34c) | 
Figs. 5- 7 are graphically ; obtained by the anmanmmmnsenes the 
_ foregoing relationships. The actual wave height at the frontal side is — 
@ 0,050~0,200 laminar, writers 
0.035~0.088 lominar , Mayer 
048 4 Binnie 7 


xv 
" 


Vi 


3 


7 
7 less than tl the calculated one, because the calculated wave 


7 profiles have sharp 
~ wave crests due to the discontinuous shock front, whereas the actual wave 


crest” is round probably because of the capillary effect and the existence of 
-vorticities at the frontal side. 


= waa hydraulic characteristics in flow, due to the formation of roll waves, 


> 


the | details of hydraulic : characteristics | in flow will be seen in the writ 


ers 

‘Based on the results derived and the obtained by the w writers 

the research works, the following 16 items are presented: 
_ 1. The author states that the roll wave is the result of interaction between 7 

the surface tension and the gravity force in a slightly disturbed ‘Soe . How- 


— 
= 


‘Sane biog as presented in the foregoing, without introducing the effect 
of surface tensions, and the results derived theoretically are satisfactorily in 
_ agreement with the experimental» ones obtained by the writers and _ Binnie. 
From this fact it may be concluded that the effect of capillary forces on the ~ 
a ’ formation and the characteristics of roll waves is quite little or rather negli- 

f gible practically, except in cases” 3 of steep slopes and of 
author concludes that slug result from instability which 7 

the transition from | laminar to turbulent flow. ‘Slug flows, 


low Reynolds number in flow just changed from laminar to com-_ 


= turbulent through transition in its flow regime. _ The fact that roll waves 


turbulent 


NUMBER, 


 pesateane: as described in the foregoing. Slug flows, therefore, . should be a 
‘plained as kinds of roll waves which result from instabilities of flow itself, at 
low Reynolds number in turbulent flow as well as roll waves in laminar flow. 
The author’ s Eq. llis questionable, whether it represents the 
for formation of roll waves and slug flows ornot. According to Eq. 11, roll 
‘waves are formed in laminar flow, when Np <2 . However, the author’s Fig. 
12 indicates that roll | waves were observed even when Ne > 2, and the writers 
also. obtained the same results as shown in the figures illustrated in the fore- 
going. Consequently, Np = (2 is not a criterion for formation of roll waves in 
author applied the criterion for the of turbulent flow and 


— 
— "282208 
FIG. $.—-RELATION B 
a 
= 


_ which was —e by the use of the Chézy formula, to laminar flow. 
_ however, Eq. 35 was deatred by putting C = constant, - author’ 5 Eq. 21 
ee 5. In regard to the effect ‘of surface tensions, ‘the author hor suggests ¢ 5 that the all 
of roll ‘waves is ay 
; 


This relation is based on ‘the idea that wave length | decreases” 
Beeston the wave height increases with time, and ultimately the condition of | 
minimum celerity Corte is attained. _ However, it is evident from Figs. 20 and 
23 that the wave velocity, the wave length, and the wave height actually in- 
“ r rease until the terminal values are reached with the lapse of time. ali alll 
is inconsistent w with the author’s ‘Suggestion. 


laminar 


4 


be. 


= 
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FIG. 6. _RE LATION BETWEEN wae PERIOD AND 


For laminar flow, the author’s Eq. 19 the plot of Np/s! 


against Rp shown in Fig. 14 is correct as eonatentes again in 7 18. 
ame manner, since Eq. 19 is written as aS a. 


— 


const. v, y, o/p) (Nw 91/9)9/10 


‘The author’ s ‘Eq. 15 is suggested to plot N,/s!/2 | against 
a Graphical expressions for Np against Rp in Fig. 12 and for Ne satel 
Nw in Fig. 13 are the indication of the laws of resistance in each regime of © 
flow, and the hydraulic characteristics of roll waves and slug flows are not 
shown. Therefore, the author’s conclusion is not derived from these figures, 
that the surface tensions and viscous effects predominate in the phenomenon ~ 7 
s roll wave formation, and at vo tag Reynolds numbers, corresponding to the de 


— 


= (36) 
— 
o3 ait _ turbulent Ly 
— 
— 
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Froude from ‘the relation of Eq. 19 because the 
‘eee of roll wave characteristics are functions of the Froude number, as | 


formula is used and its roughness coefficient is independent o on n the 
_ viscosity. Actually, however, the viscous effect may exist more or | less ‘ift the < 
channel bed condition isa a hydraulically 


If Np and 


it is found after. computati n using tee 19 and 37 that the 
Value of 6 ‘should be 1/4 for flow instead of 0. in the author’s- Eq. 


‘FIG. LATION BETWEEN DID DIMENSION LESS WAVE HEIGHT AND > FROUDE 


‘ 9. tn obtaining the minimum slope 0.029 pon the f formation - roll ag 


was ‘This c eriterton i is corresponding 


DISCUSSION 


= that the latter is the author’s idea for the roll v wave 
(2. In deriving the minimum slope 0.01 for the formation of slug flows, | 
5 Eq. 39 was also applied. At the Reynolds n number 1,200, , which the author used 
_ asa lower limit of Reynolds numbers, the flow regime will be the transition © 
from laminar to turbulent flow. Therefore, it is doubtful that Eq. 39, which © 
_ was derived by using the law | of resistance to laminar flow, c an be applied to 
the formation of slug flows. 
11. The writers are much interested in Fig. 19 which shows the location of 


initiation of roll waves and slug flows. This figure indicates the fact that roll 

waves and slug flows cannot be formed in the transition from laminar to tur- 

-bulent flow, otherwise very long distances are needed for their r development. 


‘This fa was al out nad writers in their Fig. 1 


= 


8. OF THE WRITERS’ AND BINNIE’S DA DATA. OF ROLL W WAVE 
FORM PROPOSED BY THE AUTHOR. 
As 12. The author suggests that the wave velocity of roll waves Vw is equal to i 
min = 0-763 ft per sec. Although the value of Vw ,-V - 


in which C 
tained from the experimental data is of the same order in magnitude 2 as sCmin 


his suggestion to the wave velocity is too intuitive and the effect of captlincy 

: 18. The writer’s and Binnie’s data of the wave velocity in laminar flow in mS 

presented in Fig. 8 in the same way as plotted in Fig. 24. It is found from ar 

figure that the author’s data are larger than the writer’s and Binnie’s, a. 


@ Reynolds numbers of 280 to 420. _ This disagreement in addition to ) that in ~ 
Fig. 4 cannot be understood by the writers. 

The ratio of the wave velocity V,, tot t 
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= of the slope S. Therefore, the slope should be introduced as a parameter in = 
j Fig. 25. It is not difficult to! find from Eq. 42 that the value of Vy /V is isalways 
greater than unity and increases with decrease in Rp if the slope is s constant. 
This trend agrees with that shown in Fig. 25.00 
-- 15. In turbulent flow, the ratio of the wave velocity (slug flows) to the sur- ; 
is expressed, from Eq. 6 if assuming a = as 


_ Although it seems possible from Eq. 43 that the value of Vw/Ve becomes less 

than unity, which depends on the value of V/Vs, , the results of his observation, 
which indicates that the surface velocity uceeda the wave velocity of slug | a 
“flows, can not be understood. If so, when the flow with waves is transferred to 

_ steady- -state we moving the coordinate system at a constant wave velocity, the 


velocity in downstream direction. This phenomenon is possible only in in 
vicinity of the ‘wave crest > ieee from ™ horizontal vortex due to a 


16. It may be rather better than the plot in Fig. 26, that the ratio of the re 
wave height to the wave length D}/L is plotted against Froude number with a i 
parameter of the slope because both Dp g cos and L tan 6g cos 

are functions of NF as presented previously. 
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SPILLWAY DESIGN FOR 

Discussion by B B. Michel and A. A, R. 


8B. MICHEL” and A.R. GAGNON,* A.M. ASCE.—This paper has brought up 


“one of the oldest, but still one of the most discussed, problems of hydraulic — 7 

engineering: the design of spillway crests. 
_ The author should be commended for having presented his interesting views 

on this problem é and for the attempt to standardize the design of spillways. —- 

as In this last regard the writers agree fully with Mr. Webster that semicir- 

cular pier noses induce small flow contractions and that a 1 on 12 taper on the 


- downstream end of gate slots | is generally one of the best means to lower cavi 
_ tation risks at the slots. | 


_ However, the writers believe that a second thought should be given before 

generalizing the process of designing: spillway cres cress for for 15% of the nominal 

An contribution in that field been 


spillway crest design. = Mr. Lemoine proposes to adopt a cavitation factor for 


spillway crests in the same manner as for the design of turbines: — 
‘This factoris 


in which py is the water v vapor pressure at the site, p ‘represents the absolute ; 

mia pressure on the spillway crest, v< 2/2 g is the velocity head at the 
point of minimum peeneure, and w denotes the specific weight of water. = hale 
_ This factor takes into account, not only ‘actual permanent negative 
pressure in the dangerous area, but also the danger of having the velocity head 

transformed into an additional negative pressure in low pressure nuclei ed- — 


4 dies, formed behind asperities of the spillway surface. . Thus, the factor should 


cover all conditions of cavitation on spillway profiles. .-- hy alin oe) 


art The design of any spillway can then be made — the > following ¢ considera- 


certain cavitation; 


am a August, 1959, , by Marvin J. Webster, 

a Prof. of Hydr. Engrg., Laval Univ. _ 


Lecturer in Fluid Mech., Laval Univ., , Quebec City. 
3 “Cavitation on Weir | Crests Fs sacl R. . Lemoine. Pre from 6th Meeting of the 
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+ cavitation probably impossible, even if all the velocity head ¥ 
_ transformed into additional negative pressure; and 


5: acceptable risk of cavitation if this working manne occur 
- rarely and if the spillway surface is very smooth. 


wil 
oe the cr factor. 7 It is indeed interesting to note that the foregoing refer- 

‘ n ences show that, for an overdesign of 75%, the cavitation factor, after a “ cri- 

tical design head,” is invariably zero while the minimum pressure becomes 

equal to the vapor pressure. It might _ be also noticed that the c classical form 
‘is not generally the one that will give the highest discharge coefficient for a = 
certain value of the cavitation factor. Model tests must be made in each case. ~ 

model results given in Mr. Webster's paper, for the Chief Joseph 
ee show a minimum pressure of 23 ft of water with no piers and -9 » 


the fact that very small differences bos ‘the geometry of each project have an 
important effect on pressure values, and also by the very high difficulty in-_ 

— volved in measuring the peak of rapidly varied pressure distribution on a 
scale model. The writers do not believe that the results for this spillway y are 
of all those of the same type. 

_ Taking the last value for the minimum pressure of -9 ft with an 30, which 
head of 54 ft it is found that the cavitation factor would then be 0.38, which is 

_ somewhat lower than could be permissible with Mr. Lemoine’s theory. © 

‘The author was fully aware of some cavitation risks and states that due to 
the low predicted frequency of the maximum discharge, the cost saved in in- 
_itial construction more than offsets the cost of infrequent repairs” — 
_ from damage due to negative pressure. 
The writers agree with him and this’ could easily be taken into account by 
choosing, in each case, anaccepted value of the cavitation factor. The writers _ 
strongly believe that this would be more logical than to accept a percemage 
head design of the value of the — shape 


4 “Untersuchugen a an 1 Uberfallen,” by O. Dillmann. . Mitteilungen des hydraulisch en 


5 “Engineering Hydraulics,” Hunter Rouse (Editor), Chapter VIII. ck 
ow 6 “Model Research on Spillway Crest,” by Rouse and Reid. Civ. Engrg., Vol. 5, 
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CONTROL ASPECTS OF ‘CAUCA VALLEY DEVELOPMENT? 


tng 


"Discussion by Julio Escobar- Fernandez 


Lope? 
ESCOBAR- FERNANDEZ. 1_Three aspects” 
Messrs. Kirpich and Ospina could be considered. First, the significance and 
importance, ‘in the development ¢ of the so-called underdeveloped countries, of - 
the establishment of regional autonomous corporations of development; 
ond, the technical difficulties that have to be faced by engineers interested in — 
water problems in such countries, where hydrological data is “invariably 
: scarce, or of doubtful value and inconsistent; and third, the technical features a 
7 ofa big hydroelectric project, fleetingly mentioned in the paper, which could — 
be developed by diverting the waters of the Cauca River to the Pacific —_ 
a a a system that could be classified as one of the important energetic possibili- 7 
ties, not only of Latin America but of the world, with the particular advantage 
of its proximity to the sea, a feature that makes it especially attractive for 
the establishment of heavy electrochemical industries. 
_ Although the three aspects might be of general interest, , due | to limitations — 
of space, only the first one will be discussed, at present. Some of the general | 
_ principles governing the policies and the organization of regional corporations — 
devoted to local development are set forth, and some emphasis is placed on 
_the importance of such bodies as helpful instruments in the execution of wider — 
development programs, usually national in scope. The timeliness of the 
seems especially relevant today (1960), “when better management for 
economic improvement constitutes a daily concern for the governments in the 
_ Colombia is a country with rich natural resources that can be developed to 
‘create an economically solid nation, thereby raising the standard of living of 
“its people. _ To reach these t two ) objectives, in time and without waste of effort, 
it is considered that a plan for the development and effective use of those re- 
££ rces is indispensable, as is the implementation of the plan. In Colombia, 
the Central Government has ; already launched such a program and the recent- 
me ly organized Planning Unit, directly responsible to the President of the Re- P/ 
* public, is a fundamental step towards the coordinated execution of the plans : 
that have been already formulated or which are under study. Asa masting ; 


primary policies of development, because such investments as in transporta- 
tion, energy production, irrigation and are or should be the 


cin 


foundation on which a solid economy must rest. 
On the other hand, the speenen 2 of public enterprise, which constitutes the 


Escobar, V 


| 
a 
1959, by Phillip Z. Kirpich and Carlos S. ‘ 
enegas y Rodriguez, Ingenieros Ltda., Bogota, Colombia. 4 
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stimulated not ois by the e facilities provided by those public investments, but © 
by the creation, through government measures, of a minimum of general con- | 


‘establisha “balance between all resources available— -economical, physical 
-human—and the many needs of the country. From this overall picture, a more > 
or less rigid schedule of priorities must be determined, in order to establish 

= and v where such resources can | be spent ‘most usefully and effectively. — = | 
— Once a a general plan for the whole country has been defined, the next step 
should be to project it on a sectorial basis, ‘to use the economist’s jargon, or 

on a regional basis, if the characteristics of the country so allow, as is the — 

F OA program of integrated ‘regional development, such as the one undertaken 


Colombia the of the Cauca River has, others, the 


he e establishment of a kind of experimental field sev-_ 
“eral techniques of development can be tested. In this way, the selected zone a 
"becomes a demonstration area , where e direct experience is is the best ‘school; a 


nation- wide program of development. reduced; andy 
administrative problems are also reduced; and 


out the country. 


When in 1954, David Lilienthal, mentioned by tl the authors of the paper under — 
- discussion, recommended to the Government of Colombia the creation of an _ 
- autonomous regional authority,2 as a first stage in its wider program of de- i, 
velopment, he based his recommendations on several general principles 
; should guide this type of body, principles that can be resumed as follows: ‘aon 
The initiative fora general program must originate in the region itself. It 
- ‘therefore indispensable that people in the area be willing to cooperate | in 
_ the plans, and that enthusiasm for innovation and improvement be general =, 
; ie Second, the approval an and ‘cooperation of the Central Government is s also 
"necessary. 4 _ The Government should be disposed to frame the program with 
- some autonomy, through the creation of an agency or public corporation de- | 
voted to that purpose. This agency should be provided with th sufficient admin-— 
istrative authority, and with well defined functions and aims. To better saad 
the status of that corporation, it may be useful to quote the words of President _ 
Roosevelt when he asked Congress to create the TVA in 1933 (quoted by Mr. 
Lilienthal, : : The authority should be, he said “a | corporation clothed with the = 


Fourth, the creation of a local economic corporation to 
mote and stimulate private enterprise in its own development programs. __ 

_ Finally, it is advisable that the selected region for an experiment of this 


sort, especially if it is the first to be undertaken in the country, should be one 7 
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Where the benefits obtained through the techniques of planning, are more or 
: _ less immediate, in order to overcome the natural tendency of people to be 
skeptical or even hostile to programs that require a long time to produce — 
All of these conditions were satisfied by the Colombian region of the Cauca 
‘Valley, and the CVC, the autonomous corporation created in 19 55, ts —— 


iz to show results in a has that may be considered a real s success. 
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Discussion by A. Chamberlain and Nicholas jas Bilonok 


1A, M. ASCE.—The are to be for 
_ their very careful work in obtaining design information on friction factors for 
large-diameter corrugated pipes. _ These data and ‘results presented by the 
ee will undoubtedly be used by design engineers for several decades, in 
_all cases where pipes of thenature tested are used. This paper represents the 
result | of a tremendous amount of effort on the part of Army Corps personnel. - 
‘The writer wishes: to extend the data presented by the authors on the basis 
of clear water, full pipe flow data obtained during 1953-56 in the Hydraulics 
Laboratory, Colorado State University, Fort Collins, Colo. These data were 
_ obtained incidental to certain basic research studies on nominally 12-in. di- 
‘ameter smooth, standard corrugated and helical corrugated pipes. 
_ Each of the types of corrugated metal pipe sections tested were 100-ft long, 
> ‘set on a zero slope. _ Both corrugated pipes had a clear inside diameter of — 
1 + 0.03 as as determined from ten trials on each pipe. A 
The water was recirculated ina a closed system. The discharge \ was meas- 
ured by means of a sharp-edged orifice which had been calibrated volumetri- 
cally while in place. Discharges were determined, for eachdischarge-control- ¥ 
 : valve setting, from the orifice calibration curve using the arithmetic average 
j of four separate readings of the differential head across the orifice. © 
ae. Hydraulic gradient was determined by means of water piezometers located 
at 10-ft intervals along the pipeline. The piezometer openings into the pipeline 
Bets 3/64-in. in diameter. For each discharge, each piezometer was read 
Sour times and the arithmetic average used in hydraulic gradient computa- 
tions. . Upstream « end (for 20 pipe diameters) and downstream end (for 10 pipe a 7 
diameters) piezometers were generally ignored, though recorded, in hydrau- 7 
lic-gradient studies. The piezometers were located on the crest of the cor- . 
ia (viewer inside the pipe) after preliminary investigations of the ef- . 


- _ The procedure was to first read the orifice manometer, then the piezome-_ 


ter bank and next the orifice manometer again, until each had been read fo four 


times during a period ‘of about 1 hr. Next the discharge valve : setting was 


changed, 15 to 20 min. allowed for the flow t to stabilize and the next set of 


C1 shows the two boundaries tested and reported in this discussion. 


standard ‘corrugated pipe was zine coated close- riveted construction. 


September, 1959, by M. J. Webster and L. R. Metcalf. 


Acting Dean, Coll. of Engrg. and d Chi, , Engrg. Research, Colorado State Univ., Fort © 
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tude of the corrugations from the n mean was 0.0208 ft. . The « crests and troughs 
were circular arcs with an included angle of approximately 53.73 deg and a 
radius of 0.0573 ft. Straight tangent sections joined the circular arcs. aes 2 
_ The helical corrugated pipe was zinc coated. It had a continuous lock seam 
joint. The coreagntens See a pitch or wave length of 0.167 ft normal to the | 
axis of the corrugation. The amplitude of the helical corrugations was 0.0185 
ft, measured from the: mean elevation. To complete the corrugations, circular 


ite 
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length ‘in inc 
Section views are to. ‘cor 
—DETAILS OF BOUNDARY FORMS. 


helix lix angle, ‘measured drawn along the outer. extremity of 
pipe parallel to the pipe axis and a line parallel to a corrugation, was 66 deg. 
ns aay Garde’s data were obtained with clear water full pipe flow after the pipe- 
lines had operated several hundred hours, transporting 0.2mm mean diameter 

sand at varying concentrations. Chamberlain’s data, designated by a cross, on — 

the Helcor pipe only, were taken just prior to Garde’ s. . The remaining data, 
all taken by Chamberlain, were recorded with clear water full pipe flows an 


Derived quantities Re, f, and V1 were either computed directly or read from 
graphs prepared for the purpose. In either case, the pipe diameter employed | 
in the analysis was 1.00 ft, slightly less than the inside clear diameter. Piuet ; 


magnitude of of the. diameter was used because (a) it was very convenient, 
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more important it was ‘felt that most ‘engineers applying the data 
automatically use 1.00 ft in all their design computations. 
; - Tables 1 and 2 summarize the full pipe flow data obtained, admittedly not bo 
very extensive, but hopefully of some use to the profession. These data are b 
_ presented in Fig. C2, for both the standard corrugated and helical corrugated, 
in a manner comparable to the 
author’s Fig. 10. They are also as 
- given in Fig. C3 in the same | way 100 FT. DIAM. 
Fig. C3 indicates that a rec- 
ommendation for new 1-ft diam- bg 
5.0 
eter ‘helical corrugated pipe, as 
= 0. 040, is appropriate. This 
figure shows that f isnota 5 
function of Reynold’s number a 
be partially explained by the tact 5 
the corrugations are 
to the direction of flow rin 3 
thepipe. 
‘The writer is unable to say = 
why Garde’ s data, previously = 
published, have such a wide | scat- 5 9.40 
ter. This is particularly true 
the helical corrugated pipe data — Sea 
since the pipeline was not in use 
during the interval between 
‘Chamberlain’ s and Garde’s used-_ 
Perhaps the most important *Chombericin 
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point to note is that the dashed a aa 
portion at “the smaller diameter 


extended by plotting in for stand- 


"Discharge in 


Ls FIG. C2. C2.—DISCHARGE VERSUS HYDRAULIC 


ard corrugated pipe the tentative- 
_ ly recommended point f = 0.12, at 
= 1.0 ft. Also to be plotted GRADIENT FOR 1-FT. DIAMETER HELCOR > 
for corrugated pipe is the CORRUGATED PIPES. 
n=0.0255 ata pipe ‘diameter of 1 ft. The data are insufficient to if 
Fig. C4 graphically | describes ‘the variation of the total differential head 
between crest and trough for the standard pipe tests 


of the author’s Fig. 17 can be 5 


ae also in Fig. C4, of the ratio of the hydraulic gradients as determined by 
piezometers at the crest(s) | and at the trough(s) of the corrugations, as viewed | 
from inside the pipeline. _ This plot substantiates that the author’s and the 
writer’s method of determining hydraulic gradient from a ‘set of piezometers, 7 
located consistently at the same relative location on the corrugations along a 

the pipeline, is sufficiently accurate for e engineering Purposes — the ; range 
_ of Reynold’s number incorporated in the tests sts. a 


— 
‘as 
— 
4 
— 
q 
— 
be Be 


March, 


TABLE 1. _—SUMMARY OF TES TEST RESULTS FLOW 


FT-DIAMETER CORRUGATED PIPE® = 


feet 
r second 


655 0.131 
7 


a Data from R.J. Garde, “Sediment Transport Through Pipes,” Thesis in partial ful-— 


fillment of the requirements for the Master of lorado Univer- 
Fort Collins, Colorado, December 1956. 
a Taisen Inch Pipes,” Dissertation in partial fulfillment of the requirements for = 
Doctor of Philosophy Degree, Colorado State University, Fort Collins, Colorado, June — 


«5.98 6.80 | 9.0 Average | 0.117 
? 
| 
‘ 
“| 
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DISCUSSION 
TABLE 2.-SUMMARY OF TEST RESULTS FULL PIPE FLOW 


j inecubic feet | 
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Data from R.J. Garde, “Sediment Transport Throu | the 
_ requirements for the Master of Science Degree, Colorado State University, Fort Collins, Colorado, 
Db Data from A.R. Chamberlain, “Effect of Boundary Form on Fine Sand Transport in Twelve-Inch _ 
—_ Pipes,” Dissertation in partial fulfillment of the requirements for the Doctor of Philosophy Degree, 
Colorado State University, Fort Collins, Colorado, June 1955.00 
© Data obtained by Chamberlain after pipeline had been used to transport 0.2 mm san 


Arithmetic average, all Chamberlain data, f= 0.0400. 
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FIG. C4. RATIO oF ‘HYDRAULIC GRADIENT on CORRUGATION HIGHS TO THAT OF 
CORRUGATION LOWS, AND DIFFERENCE IN HEAD BETWEEN HIGHS AND LOWS 


= presentation of data from certain carefully conducted experiments on full _ 
size, , large commented metal pipes, to determine friction factor “f” abso- 
lute ‘Toughness n.” The accuracy depends on a correct evaluation of the 
well known coefficient of absolute roughness in the Manning formula. aa 
ae depends primarily < on the surface c condition of the perimeter of a pipe or open- — 

Because of the co complex nature of the friction factor various 
have been made to group all the roughness dimensions together into a single 

‘coefficient of absolute roughness “ n” for simplicity. 

a. The writer is s particularly interested in this paper as he had independently — 

investigated the relationship between friction factor “f” and absolute rough- 


Civ. Engr., Engrg. Div., of Engrs., U. Ss. Army, Philadelphia, Pa. 
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values of Q, S, H Hy, L and D have been measured in and 
friction fac factor “f” has been found be be given | by 


a 


a 
R2/3 g1/2 


3 may be written from the friction factor by the 


-@ 


n” ranging from 0.010 to 0.025. 


‘The values of “f” computed from Eq. 6 a are in those 
= by experiments for —- = ft, 7 ft, and 7 ft in 1 diameter, and also with © 


— 
given 


The exact form, and position Reynold’ s number for 


= practical value of ‘the aye number is that it indicates the degree _ 
 e turbulence in flowing water. Fora given size of circular pipe the velocity — 


is the major variable the Reynolds number increases | as the velocity of 


— 
— 
if 
— 
a review experiments data presented in the paper the writer has com 
ness coeffi- a 
— 
— 
a 


1960 


velocity of flow ‘The water changes have no significance. 
Investigations | of friction losses_ for Part-Full flow, in Reference 3, 
found that ' “For uniform oube ritical flow, , the Manning coefficient did not indi- — 


SS: 
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i. 
1 
0 
0. 
9. 
9 
8 
8 
8 
8 


cate: a systematic variation wn: Reynolds or pipe size the 


a — can be recognized that several factors might have affected the accuracy © 


J 0.010 0.011 0.012 0.013 | 0.015 0.021 | 0.024 |o028 
6 23.3 39.4 | 75.5 | 134.2 | 145.6 | 
30.5 | 35.7 | 47.6 | 68.5 121.8 | 132. 
fi 10 19.7 | 23.8 | 28.4 | 33.3 | 443 | 63.8 1135 |131. ii 
1 106.6 | 115.6 
98.9 | 1074 
93.4 
78.3 | 85.0 
36 74.0 | 80.0 

48 67.0 | 73.0 

72 «| 58.5 | 63.5 |, 
— | 53.3 | 580 
49.4 | 53.5 
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Discussion by ‘Steponas Kolupaila, Israel 


and William C. Peterson 


STEPONAS KOLUPAILA. | 
problem of the correction factor, known as a ‘Coriolis’ - coefficient, is still 
: active, and still not entirely clear. _ The purpose of this factor is to ‘correct 
- the ‘velocity head in the Bernoulli equation — for different velocities across the 


‘a The basic expression given under 10, p. 72 is correct, euvegt for an ob- | 
ious typographic omission: AQ. is to be understood instead of Q in the nu- 
merator. This formulacan give an accurate solution when a stream is divided 
into a large number of small areas with measured velocities. This work is 
‘more reasonably performed > by a graphical integration with a planimeter, : 
those “astronomical” ‘numbers. Application of this formula for 
several parts of a cross section with assumed coefficients is a half measure. | 
‘The value obtained by the authors (1.19) seems to be too low for anopen —. 


nel with overflooded banks. ake - The river gagings by two point | 


cient. Two integrations are necessary in an open channel, along the verticals, 
and across the channel, and more points must be taken on every vertical. —_ 
ISRAEL STEINBERG, G,3 F. ASCE.—The sample computations presented 
in the article, and which are patterned after a procedure given in the Engi- 
= Manual of the Corps of Engineers, still contain some of the e objection- 
able, relatively large, figures required to obtain a relatively small value of 
velocity head. For this reason the writer has adopted a number of -qetaiens il 
_ Procedu res in order to further reduce the necessity of dealing with large and 
unwieldly values. One of the es utilizes a weighted area curve for wall 
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Ja 2 “Open-Channel Hydraulics,” by V.T. Chow. . McGraw- Hill, 1959, p. 28. 
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in which Aw is the weighted area, Ac and Aj are respectively, the channel and © 
overbank areas, and and Co are, Feapectively, the channel and 


in which C is the total conveyance factor for the section (= =Co + Co), as Lis 


_ Eq. 1 was developed for flow lineproblems involving separation of the sec-_ 
‘tion into two parts; namely, channel and overbank. The equation can be readily 


expanded to include any desired of components. Taking the reciprocal 


a ma. 1 and squaring wi 


_ Eq. 3 applies, also, number | of ‘components a can be be written in the 


Eq. 4 can be further modified in more convenient forms for computing — 
- either (a) the weighted velocity head, hy, or (b) the alpha coefficient, a, in 
of total section area. To obtain by, both sides of are e multiplied 


of Eq. 5 velocity head is illustrated 


P in Table 1 using the author’s values for section 2. The values in Cols. 7 and 8 
were obtained with of the slide rule. le. 
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In Table 2 only Cols. 7, 8 and et are a 


the assumed value of Q = 10 
hy = 1.424 | = a 
WILLIAM PETERSON, 4M. ASCE -—Computation of an approximation of 
the velocity head correction factor, a, by the fraction = _i.... 7 
if is essentially a simple slide-rule operation, quickly accomplished when: oo 


becomes familiar with the “feel” of the magnitude of the values involved. 
Techniques, which are largely personal, abbreviate the consequent “ as- 


figures of the computation. 


— "0.006 | (9.000 0,003 
Computations are the same as forthe 0. 
| 
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the of the author’ S paper the expression for friction loss 


4 


‘The a intnahes neglect to state, however, that the quantity in parentheses must 
; 4 always be positive. For an an expanding reach, hy, > hy. vo» the ,, the quantity in in — 


‘ However, for a contracting oe hy. < hy, ,, and the quantity» in iaiainen 


hg = 
Expressing an of a in terms of section K ‘and A, 
g essential toa direct method of determining the discharge of a stream m 
2 slope-z the two-section reach 1-2 in Fig. D1. 


= 
— 
4 ga 
4 a 
— 
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"DISCUSSION 


From continulty 


K; 812 = Ky S 
we 


he = (by, - hyo) m) ae 
f + (by, - by, 


and = he = 


— 
— 
— 
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“Modification o: of the area relationship by values determines whether 
reach is hydraulically contracting or expanding; that is, whether the weighted 
velocity I head is increasing or decreasing in the downstream direction. If 
‘ _ quantity in brackets is positive, the reach is contracting; if negative, the — 


‘Using data from the author’: s example, in the formula for a contracting, 


was ‘computed as 10,700 cfs. 
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and, for an expanding reach — 


Re 


_FRED W. BLAISDELL, 1 F. ASCE.—Mr. Moore has presented interesting 


- Oklahoma and Arkansas during the floods of April, May and June, 1957. — Ac- : 
cording to Mr. Moore, the performance of the structures, designed for a 50- 
yr storm and subjected to greater than 100- yr storms, 
_ tations. . uh The paper is devoted primarily to the performance of the emer 
- gency spillways. . This is well because it is noted that “. . . the use of vege- : 
tated earth spillways for flood prevention structures . ee a material diver- 
gence from the accepted principles employed by th the engineering profession in 
the principal spillway, Mr. Moore notes that it is used to di 
_ charge water “...ata rate designed to prevent flooding of the stream chan- 4 
; nel below.” A section through a principal spillway is shown in Fig. 3. ‘It con- 
g sists of a drop inlet, a barrel through the dam and a cantilevered outlet. _ The be 
fact that the performance of the principal spillways is not mentioned is prob- = 
ably due to the fact that they performed as designed. The writer has been con- ‘“ 
cerned with the hydraulics of the principal closed conduit spillways since : 
1941, and naturally was extremely interested in seeing how they came through: 
these severe tests. In September, 1958, he visited a number of the structures. 
The visits confirmed the supposition that the structures did their intended job 
with 
The practice of omitting stilling basins and simply cantilevering the pipe 
beyond the toe of the fill, as shown in Fig. 3, is a practice that is also “a ma- . 
terial divergence from accepted principles.” It was the performance of these 
_ cantilevered pipe outlets that the writer was particularly interested in evalu- -_ 
i ating. Compared to the size of the pipe, the storage between the crests of the 
— spillways and the emergency spillways is large. At many of the 
structures, this storage pool was partially drawn down and refilled several ee 
__ times as a result of successive rains and many of the principal spillways ran 
capacity for long periods of time. This gave maximum for the 
‘development of the scour holes anticipated toformatthe pipeexit, 
The performance of the cantilevered outlets was much better ton . 
writer had expected. Typical scour holes are shown in Figs. and 
_ Fig. B1 is the outlet at Honey Creek Site No. 8E. There is a 17-in. concrete ia 


pipe through the fill and the cantilevered pipe is 18-in. corrogated metal. The — 
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outlet ran for six weeks at full lnidibianti the 12- “in. . outlet pipe shown in 1 Fig. —— 
B2 is for the Honey Creek Site No. 14. The depth ‘of the pool in the scour » hole = 
‘is about 4 ft. It can be ‘seen in Figs. B1 and B2 that the scour holes are small 
in size. At both sites, as well as all other sites seen, the scour hole termi- 7 
nated close to the end of the pipe. The bent shown in Fig. is intended to 
_ port the pipe if the scour hole progresses upstream. The experience of the 
Texas structures indicates that this feature may not ‘be necessary under the > 
conditions existing there. However, it is not meant to imply that bents would 
u- _ The conclusion to be drawn from the Texas experiences is that a canti- 
levered. outlet may be used in place of a stilling basin 


The v writer has been lo looking at cantilevered outlets throughout the Midwest oy 


for a number of pense with the expectation of finding some structures in which 


cures in which the scour holes looked as oes they might eventually _ 
_ develop to a damaging size but which hadnot yet progressed to that point. One ~ 
_ of these outlets is shown in Fig. B3. This is a 42-in. cantilevered outlet built 
4 1950. It carried a designed flow in 1952. _ This single flow formed the wide 
scour hole shown in the picture taken in 1954. Comparison with a picture — 
taken in 1952, , shows no visible increase in scour hole dimensions, — Ps 
because additional capacity flows have not passed the structure. a Both pic- 
- tures show ‘that the scour under the end of the pipe is insufficient to pene ; 
7 it in any way. In fact, the scour has not undermined the pipe as far as the bent | 
provided to support the pipe | over the scour hole. In spite of the rather wide 


_ scour hole, the outlet performance must be classed as satisfactory. Whether 
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bag: writer feels unsure al the performance of cantilevered outlets in spite 
of the fact that he has been unable to discover field evidence to back up his — 
doubts. — a He feels there must be limitations to the size of cantilevered —— 
that can be used with soils of different erodibilities in order to prevent the | 
scour hole from endangering the dam. Here isa fertile field for research! 
will the scour holes develop laterally, rather than in a dangerous upstream 
direction? Will upstream development of the scour hole take place : after the | 
lateral development has reached acertain extent? How long will it take to de- : 
velop the scour hole? What is the effect of the bed material? What will be its 
ultimate size? Is this an economical self-formed stilling basin? These are 

— 
_ Evidence collected by the writer indicates that this type of outlet structure — 

; ni a definite usefulness. | Lack of difficulty with it is such that there is no | 
pressure from designers to conduct a thorough investigation into the charac- 
teristics and limits of applicability of the cantilevered outlet. However, the 
writer feels that research on the cantilevered outlet would be of of major value 


“a -__ or not future capacity flows will increase the size of the scour hole to dan aa 
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VORTEX CHAMBER AS AN 
-FLOW- CONTROL DEVICE* 
Discussion by Michael Amein 


_ MICHAEL AMEIN IM. ASCE.—The formation of vortices in flow, in “ 
structures, mas been found of three principal 


Ca cost for the structure. The third feature is the possibility of damage 
_ from the erosive forces of the vortex motion. Such forces have been known to 
{ 7 ‘cause damage to inlet drop spillway structures at earth dams. * The authors» 


a case in which the presence a is not accompanied 


rortex in chamber. the applications | by on ero- 
_ sive forces are either insignificant or irrelevant, and the reduction in = 


-Pacity would be used to advantage. 


type vortex. described by the authors usually ‘results from flow 


derived by the application of Bernoulli’s along and by 

equating the radial pressure differential to the centrifugal force. In = 

analysis the effects of the viscous forces are neglected. _ However, since in 2 
real fluid viscous forces are” ‘present, the: free vortex theory ‘not 


orbits the core and are discharged through the orifice. ‘Thus, the par- 
ticles “constituting | the vortex are constantly being renewed. It is the pre- 
| tabianes of the tangential velocity in bees vortex chamber that gives thi 
4 
of free vortex in this case is probably correct for all practical purposes. 
The authors have presented the results of their experiments in terms a 
_ coefficient of discharge and adimensionless parameter called the vortex num- > _ 
_ ber. This parameter contains the prime factors inthe vortex motion. The 
e writer would favor the publication of the observed data in conjunction with the »” ; 


4 results. A better understanding of the physical the of 
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the head-discharge relationship, resulting in an erratic hydraulic behavior. 
— 
— 
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a 
mi _ tion as enunciated by Helmholz, the free vortex consists of the same fluid | nem 
Fe 


and the quality of experimentation could ould be the 
a observed data. Although the presentation of data in dimensionless parameters 
is a commendable procedure following the widespread acceptance of dimen- 4 
‘sional analysis, yet, this type of presentation, unless Supplemented by ob- 
served data, freezes the data to a single interpretation. Furthermore, the _ 
| data would not be accessible to. support or refute other approaches that might 
¥ 


‘The reduction in the discharge capacity of the orifice with the formation of — 


Pat the: vortex can be said to be due to the conversion of part of the static head to 
tte tangential kinetic energy in the vortex. In the vortex chamber, the intro- — 
, duction of flow at a tangential direction stimulates the formation of a vortex — 
5 and tends to make the vortex a stable feature of the flow. _ The author’s ex- 
s periments indicate that the coefficient of discharge C is a unique function of 
_ the vortex number, from which the stability of the head discharge caceuadl , 
— . the following analysis, three new symbols one wah in addition | to the 
; oS described by the authors: b = width of either of two inlet channels; 


* W = velocity in the inlet channel at the entrance to the chamber; and 6 = _ 


= 


The velocity in the inlet channel is so small that the depression of = 
er surface due to the — head is neglected. Then inflow et 


‘But, from the definttion o of vortex number 


| 


By substituting t this v of Ww in 


= 


3 
“4 
of 
— 
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. Eq. 7 indicates that the head H ina given vortex chamber can can determine | the 


vortex number and the the case in which b = 6 in., B=42in., 


if H is known, then Cc and can be readily calculated. Such a calculation is 
Tt is interesting to. compare the 
static head required to provide a given 
discharge when there is vorticity in TABLE 1. —COMPUTATION OF 
the flow to the static head required DISCHARGE 
_ provide the same discharge when the ees 
vorticity is eliminated. For the latter cubic 
case, a constant C equal to 0.686 is feet 
assumed. The static head required to rr 
provide a discharge Q be ‘deter- 43800 | 0.15 


0 47 0.18 
mined from 0 0.21 
—60.51 0. 24 


TABLE 2 


: 
4 
— 
— 

results of the computat -as the head on the orifice in 
The en from Table 2, that as the hea the vortex cham 
It is seen from to vortex motio | he 
: 
— E represents the pe 
— — 


may not be sneta as a flow control device, that would ‘maintain a constant or z= 
a slowly rising discharge under a rising head. In fact, as has been experi- 
enced in drop-inlet spillway structures, a ‘ealichenlly high head will drown 
out the vortex n motion and re-establish the normal orifice flow relationships. x 
Many cases could exist in which the vortex chamber would perform no better — 
than an orifice of smaller area. It appears that it is the range of heads and | 
discharges within which flow takes place that determines whether the pres- 5 
ence of a vortex in the flow can be permitted and whther the vortex | can be 
made to perform a useful function. _ Further investigation might indicate that — 
- even in hydraulic structures it may be more feasible, in some cases, to — 
_ stabilize the vortex than to make attempts towards its elimination. Research 
on) vortex motion would not only open new fields for its —— ee but would 


also add to our basic understanding of the mechanics of — 
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